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Characterization of Price Trajectories, Cash Flows and Operational 
Costs in Futures Markets with Monte Carlo Simulation 



Adriano Azevedo-Filho 
Elisson de Andrade 

ABSTRACT 

This research develops procedures to characterize and simulate - using Monte Carlo 
methods - price trajectories and cash flows, allowing the estimation of operational costs 
in futures markets. First, a brief literature review analyses models and ideas about the 
evolution of prices in markets. This review organizes known facts about the 
characterization of price trajectories and pointers to methods used to provide a statistical 
representation for futures prices models. Second, Monte Carlo simulation methods are 
reviewed. Algorithms describe the overall simulation process of price trajectories and 
functions of these trajectories. Remarks are presented on the estimation of moments, 
confidence intervals and probability distributions. Third, definitions and algorithms 
formalize cash flow evolution (day by day) and costs in futures markets operations, based 
on price trajectories, exchange rules concerning margins, financial cash flows and taxes. 
Fourth, operational cost components in futures markets are estimated and quantified in 
(a) a case study concerning the BM&F corn market and (b) with a theoretical model. The 
estimation considered observed and simulated price trajectories, derived from ARIMA- 
GARCH models, estimated from the observed price trajectories. Operational costs are 
shown to be large and uncertain when operation begins, being strongly dependent on the 
price trajectory, degree of dependence on borrowed money to manage the cash flow and 
exposure to taxes. Because these costs apply to both the long and the short positions of 
the same operation, they affect twice the transaction costs in a single operation. A cost 
risk concept is introduced and argued as being an additional obstacle to futures markets 
operations, not found in the literature. In an operation for a single investor (not a firm) 
and with minimal assumptions, the mean percent cost for a long or a short position, was 
estimated to range from 2.5 to 4.5 % of the contract initial price, for an operation lasting 
120 days, when annualized volatility is 25%. Mean operational costs (as well as 
uncertainty) tend to increase with larger volatilities and longer operations. Models with 
more realistic assumptions lead to larger mean cost estimates. Agents who need to 
borrow money to manage their cash flow will have larger costs than agents who are able 
to finance operations with their own capital. High interest rates, wide spreads between 
borrowing and investment rates, and high taxes applied to operations in futures markets 
are argued as major causes for high and uncertain operational costs. In addition, because 
there are differences in the cost of money and exposure to taxes for different classes of 
agents, the ones with smaller costs have more incentive to use futures markets. Agents 
with large costs can be excluded from these markets. These effects are argued to have 
economic implications. 
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Caracterizac,ao de Trajetorias de Prec,os, Fluxos de Caixa e Custos 
Operacionais em Mercados Futuros atraves de Simulac,ao Monte Carlo 

Adriano Azevedo Filho 
Eli s son de Andrade 

Sumario Executivo 

1 . O trabalho apresenta procedimentos metodologicos para caracterizacao e simulacao de 
trajetorias de precos e fluxos de caixa, possibilitando a estimativa de custos operacionais 
em mercados futuros. A apresentacao foi organizada em seis secoes principals. A 
primeira secao contem a introducao do trabalho. A segunda secao realiza uma breve 
analise da literatura e de resultados existentes sobre a dinamica dos precos em mercados, 
com o intuito de embasar a definicao dos modelos utilizados no trabalho. A terceira secao 
discute princfpios utilizados na simulacao Monte Carlo, apresentando modelos 
simplificados e algoritmos de simulacao considerados nas etapas posteriores da pesquisa. 
A quarta secao formaliza uma definicao do conceito de custo operacional e de seus 
componentes, em uma operacao em mercados futuros. Essa definicao e obtida a partir da 
caracterizacao simplificada do fluxo de caixa, que e dependente: da trajetoria dos precos 
observada, de encargos financeiros, de tributos e de taxas operacionais. A quinta secao 
apresenta resultados da aplicacao dos procedimentos desenvolvidos nas se§oes anteriores 
em estudo de caso e em modelo simplificado, visando caracterizar o custo das opera§oes. 
Finalmente, a sexta se§ao apresenta considera§6es finais. 

2. Uma dificuldade existente para a estimativa de custos de uma opera§ao em mercados 
futuros e o fato de que alguns de seus componentes - principalmente tributos e encargos 
financeiros - dependem fortemente da trajetoria de pre§os que sera observada no 
mercado, nao sendo conhecidos com certeza no infcio desta opera§ao. Outro problema e a 
dependencia dos custos operacionais a situa§ao do agente economico quanto ao seu grau 
de necessidade de capital externo e a natureza de sua exposi§ao a tributa§ao. Para agentes 
mais expostos a tributos e/ou descapitalizados, dependentes da capta§ao de recursos para 
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fazer frente a margens, ajustes, tributos, e gerenciamento de contratos de um modo geral, 
os custos tendem a ser mais elevados que para aqueles que dispoem de ampla capacidade 
propria de recursos financeiros e/ou sao menos tributados. A incidencia de tributacao 
pode diferir entre agentes economicos: pessoas ffsicas e pessoas jurfdicas possuem 
distintos regimes de tributa§ao - alguns tipos de pessoa juridica podem ser isentos de 
certos tributos (ex. cooperativas). 

3. Os custos operacionais estimados neste trabalho referem-se ao caso em que o agente e 
uma pessoa fisica, sendo obtidos a partir de procedimentos que consideram a evolucao do 
fluxo de caixa da operacao, condicionados por trajetorias de precos: (a) observadas em 
contratos futuros; (b) simuladas a partir de modelos estatfsticos (ARIMA-GARCH) que 
capturam aspectos relevantes de trajetorias observadas no mercado; e, (c) simuladas 
atraves de modelos que so dependem da volatilidade e do numero de dias da operacao. 
Deve-se destacar que os custos para outros agentes (pessoa juridica tributada com base no 
lucro real ou no lucro presumido, instituicoes isentas de certos tributos etc.) podem ser 
diferentes dos custos operacionais apresentados neste trabalho, sendo assunto de 
pesquisas futuras. 

4. Os resultados obtidos sugerem que os custos de operacoes em mercados futuros podem 
ser expressivos. Para uma pessoa fisica, numa operacao de 60 dias, sob pressuposicoes 
minimalistas, os custos operacionais medios atingiram valores entre 1,5 e 2% do valor 
inicial do contrato (Figura 5.8, pag. 75), quando a volatilidade anualizada foi de 20%. Os 
valores apresentados dependem da situacao da pessoa fisica com relacao a necessidade de 
capta§ao de recursos - o menor valor refere-se ao caso de agentes com ampla capacidade 
propria de financiamento da operacao, chamados no trabalho de agentes capitalizados ; o 
valor maior corresponde ao caso de agentes que necessitam captar recursos no mercado 
para financiar suas operacoes, ou seja, agentes descapitalizados (essas consideracoes 
valem para os proximos resultados). Para uma operacao com 120 dias, esse custo 
operacional medio variou entre 2,5 e 3,5%. A parcela correspondente ao IRPF, nesse 
caso, foi aproximadamente 1,5% do valor inicial do contrato, para agentes capitalizados 
ou descapitalizados. Quando a volatilidade considerada foi de 30% os custos medios 
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observados estiveram entre 2,5 e 3% e entre 3 e 4,5%, para 60 e 120 dias, 
respectivamente. Os resultados apresentados assumem taxas de juros nas aplicacoes e nas 
captacoes, respectivamente, de 1,2% e 3,75% ao mes. Em situacoes praticas, as taxas de 
captacao podem ser muito superiores a esse valor, nas condicoes atuais do mercado, o 
que levaria a custos ainda mais elevados. Do lado das aplicacoes, assume-se que os 
saldos positivos de ajustes podem ser investidos a taxa de aplicacao considerada (algo 
facilitado pelo FIF, um fundo que remunera os saldos positivos dos ajustes diarios). Os 
custos apresentados foram obtidos a partir de modelos simplificados que foram muito 
uteis para caracterizar a natureza do impacto da volatilidade e tempo (em dias) da 
operacao, nos custos da operacao. Esses modelos tenderam, de um modo geral, a estimar 
custos inferiores aos custos obtidos por modelos mais complexos e mais aderentes a 
realidade, tambem examinados na pesquisa. 

5. E importante perceber que os custos apresentados no paragrafo anterior sao aplicaveis 
as duas contrapartes de uma mesma operacao, tanto na posicao vendida como na posicao 
comprada. Isso significa que a participacao do custo operacional em uma transacao sera a 
composicao dos valores medios obtidos para as duas contrapartes. E possfvel ainda que 
numa operacao os custos operacionais para a posicao vendida e a posicao comprada 
sejam diferentes (assimetricos), em funcao da forma com que os encargos financeiros e a 
tributacao sao aplicados. Numa posicao de venda, uma trajetoria de precos ascendente 
pode motivar o pagamento de um elevado volume de encargos financeiros, pelo uso de 
recursos externos para deposito dos ajustes diarios. Para o contrato de compra, por outro 
lado, pode haver um expressivo recebimento de juros decorrentes da aplicacao dos 
ajustes positivos recebidos, mas os desembolsos com tributacao podem ser elevados. Um 
exemplo que destaca esses efeitos, analisado neste trabalho, e ilustrado pelo contrato do 
milho com vencimento em JAN03, onde o preco passou, em 120 dias, de R$ 16,20 para 
R$ 27,00, por saca de 60 kg. Se uma operacao fosse realizada nesse periodo, para um 
contrato de venda, os custos associados a encargos financeiros chegariam a 9,85% do 
valor inicial do contrato, para um agente descapitalizado, e 3,12% para um agente 
capitalizado (Tabelas 5.3 e 5.4, pag. 60). O IRPF pago seria nulo nesse caso e os custos 
operacionais totais seriam de 4,18% e 12,2%, para os agentes capitalizados e 
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descapitalizados, respectivamente. Nessa mesma situacao, para um contrato de compra, 
teria havido um recebimento de juros equivalente a 3,57%, e um pagamento de IRPF 
equivalente a 15,06%, com um custo operacional de 13,6% do valor inicial do contrato. 
Se os agentes mantivessem suas posicoes pelos 120 dias, teriam sido esses os custos 
observados, que foram muito diferentes (assimetricos) para posicoes vendida e comprada. 
Se um agente capitalizado, na posicao comprada, fosse isento de IRPF, teria terminado 
sua operacao com ganhos lfquidos associados aos juros recebidos. As diferencas podem 
ser acentuadas se as operacoes envolverem agentes com diferentes graus de capitalizacao 
e exposicao a tributacao. 

6. Um aspecto importante observado foi a expressiva incerteza que existe no infcio da 
operacao quanto ao valor dos custos operacionais ao seu final. Numa operacao com 60 
dias, considerando uma volatilidade anualizada de 25% e agentes descapitalizados, o 
custo total medio foi proximo de 2,5%, podendo variar, de acordo com a trajetoria, entre 
0 e cerca de 7% (Figura 5.9, pag. 76). Os resultados mostraram que a incerteza pode 
crescer quando aumenta o numero de dias da operacao e a volatilidade. Depende tambem 
da natureza do processo estatfstico que caracteriza a dinamica dos precos (para uma 
mesma volatilidade e numero de dias). 

7. A incerteza com relacao ao custo em uma operacao em mercados futuros pode ser 
caracterizada como um risco de custo , um conceito que acreditamos ser novo e 
particularmente importante para as operacoes no mercado futuro brasileiro. Esse risco de 
custo , ao lado do risco de base, pode dificultar operacoes, pelo impedimento de um hedge 
que possa efetivamente fixar o preco (e margens em operacoes de arbitragem). Esse 
problema seria menos serio se o custo fosse fixo, conhecido exatamente no infcio da 
operacao, algo que pode nao ocorrer. 

8. Ha algumas implicacoes economicas importantes decorrentes de custos expressivos, 
incertos e assimetricos entre agentes. Custos elevados e incertos tendem a reduzir o 
beneffcio de uma operacao de fixacao de precos via mercados futuros, inibindo 
negocia§oes e favorecendo a fixacao de pre§os por outros mecanismos que tambem 
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podem garantir precos futuros e eventualmente adiantar recursos: CPRs, fixa§ao de 
precos junto a traders e processadores etc. Os precos fixados por esses outros 
mecanismos, contudo, tendem a incorporar desagios que podem ser, em parte, motivados 
pelos custos e incertezas existentes na fixacao de precos pelo mercado futuro. Agentes 
economicos mais capitalizados e/ou menos expostos a tributos, podem ser mais 
estimulados a utilizar o mercado futuro que agentes descapitalizados e/ou mais expostos a 
tributacao. Essa situa§ao e de uma certa forma perversa, pois mostra que ha barreiras de 
acesso ao mercado futuro (via custos de operacao mais elevados) para muitos agentes que 
potencialmente poderiam se beneficiar dessas operacoes (ex. produtores rurais do tipo 
pessoa ffsica). Para trabalhos que estimam a razao otima de hedge as implica§oes podem 
ser grandes, em decorrencia de que em geral assumem custos baixos e fixos nas 
operacoes. Um entendimento mais completo das implica§6es economicas desses 
resultados, contudo, exigira um esforco de pesquisa que foge ao escopo deste trabalho. 

9. Deve-se destacar que, na literatura internacional, os custos operacionais em mercados 
futuros tendem a ser assumidos como sendo de pequena magnitude quando comparados a 
outros custos e, por essa razao, usualmente desconsiderados. Essa pressuposicao, apesar 
de muito utilizada (explicitamente ou implicitamente) na literatura brasileira sobre 
mercados futuros, parece ser pouco realista. Cabe discutir aqui algumas razoes que 
motivam esses custos elevados e incertos, em magnitude potencialmente muito superior a 
observada em outros pafses. 

10. Uma primeira razao deve-se a: (a) taxas de juros elevadas praticadas no Brasil, em 
geral mais elevadas das taxas praticadas em outros pafses; e, (b) diferenca significativa 
entre as taxas para aplica§ao e capta§ao. Esses dois fatos encarecem as operacoes, 
especialmente para os agentes economicos descapitalizados, que precisam captar recursos 
para manutencao de suas operacoes de fixa§ao de precos. As taxas de juros afetam os 
componentes do custo operacional associados a margem de garantia e do fluxo de caixa 
em geral. Os custos associados a margem podem, porem, ser conhecidos com relativa 
certeza no infcio de uma operacao, o que nao acontece com os custos dependentes de 
ajustes diarios, que sao incertos. As distor§oes existentes nas taxas de juros acabam tendo 
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como conseqiiencia perversa o afastamento dos agentes economicos mais 
descapitalizados. Esses agentes, em muitos casos, sao os que mais podem se beneficiar 
das operacoes nesses mercados. 

11. Uma segunda razao deve-se a tributa§ao. Em operacoes envolvendo pessoas fisicas, 
analisadas no trabalho, o IRPF e elevado e produz distorcoes. Da forma que e definido, o 
imposto acaba motivando uma parte consideravel do custo e da incerteza em uma 
operacao, favorecendo com isso a utiliza§ao de outras formas de fixa§ao de precos que 
sao menos tributadas (ou mesmo nao tributadas). Como os tributos (IRPF e CPMF, no 
caso de pessoa fisica) dependem da trajetaria seguida pelos precos, contribuem 
significativamente para o componente incerto do custo da operacao. Nao e clara a 
justificativa economica para a tributa§ao, tanto do lado da magnitude das alfquotas 
definidas quanta do lado das regras para calculo dos tributos, que podem inibir operacoes 
que sao importantes para reduzir riscos em muitas situacoes. 

12. Os procedimentos desenvolvidos podem ser estendidos para a analise de outras 
situacoes envolvendo agentes do tipo pessoa juridica, com diferentes regimes de 
tributa§ao, cooperativas ou institui§oes financeiras, que em alguns casos sao isentas de 
certos tributos. Essa extensao da analise torna-se interessante dado que sao operadores 
tfpicos nos mercados futuros. Os custos operacionais para esses agentes podem ser 
diferentes dos obtidos neste trabalho, em fun§ao de outros mecanismos de tributa§ao. A 
extensao dos procedimentos para esses e outros casos, entretanto, e assunto para futuros 
trabalhos. 

13. Finalmente, um merito importante deste trabalho e sugerir que o custo de opera§6es 
no mercado futuro, no Brasil, pode ser muito superior a valores praticos freqiientemente 
considerados (ex. 0,64%), que tendem a levar em considera§ao, simplesmente, as taxas 
operacionais (TOB, registro e emolumentos). A pesquisa apresenta estimativas desses 
custos operacionais para pessoas fisicas, argumentando que podem ser elevados. Mostra 
tambem que esses custos sao incertos no infcio da opera§ao, caracterizando um risco de 
custo , que juntamente como o risco de base e outros custos de transa§ao, podem 
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desestimular operacoes em mercados futuros. Ao considerarmos que as duas contrapartes 
da operacao estao sendo oneradas, o custo transacional associado a operacao (podem 
existir outros) e a composicao dos custos operacionais de cada contraparte, que podem 
ser assimetricos. Esforcos na direcao da mitigacao desses custos e de sua incerteza, 
podem favorecer um maior interesse pela utilizacao dos mercados futuros como 
mecanismo de fixacao de precos futuros. Iniciativas como o FIF, o uso do Funcafe para 
financiamento de margens e ajustes (no caso do cafe) e o Fundo Futuro Agropecuario 
(Banco do Brasil e BM&F), podem ser importantes para a reducao do custo operacional. 
Com relacao a tributacao incidente em mercados futuros, e diffcil entender sua 
justificativa economica nos moldes atuais. Tal tributa§ao tende a introduzir distor§6es 
desnecessarias e inibir opera5des nos mercados futuros importantes dentro do 
gerenciamento de riscos associados a varia§oes nao antecipadas de pre§os, ou mesmo 
cambio. 
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Caracteriza^ao de Trajetorias de Prec,os, Fluxos de Caixa e Custos 
Operacionais em Mercados Futuros atraves da Simulac,ao Monte Carlo 

Adriano Azevedo Filho 
Elisson de Andrade 

1. Introdu^ao 

Este estudo discute a caracterizacao e simula§ao de series de precos, fluxos de caixa e 
custos operacionais dentro de operacoes em mercados futuros. Os procedimentos 
desenvolvidos sao aplicados a estimativa de certos custos existentes nessas operacoes 
associados a tributos, taxas operacionais, margem de garantia e ajustes diarios, 
obtidas a partir de trajetorias de precos compatfveis com parametros estatfsticos que 
caracterizam a dinamica do mercado. As tecnicas apresentadas podem facilitar a 
quantifica§ao de custos dependentes da evolucao do fluxo de caixa em operacoes 
financeiras, que envolve regras especfficas para pagamento de tributos e encargos 
financeiros, facilitando a analise de operacoes de hedge e arbitragem. Um estudo de 
caso e um modelo teorico, representativo de operacoes no mercado, ilustram a 
aplica§ao da metodologia apresentada. 

Uma preocupa§ao existente e a simulacao do comportamento dos precos de uma 
forma consistente com fatos estilizados conhecidos na literatura, os quais sao 
examinados. O estudo compara o desempenho da simula§ao realizada por modelos 
estatfsticos alternatives, que consideram diferentes graus de complexidade e 
aderencia a realidade. Essa compara§ao tern o objetivo de avaliar a sensibilidade dos 
resultados aos modelos utilizados, verificando se os modelos mais simples podem 
solucionar satisfatoriamente o problema em questao. 

O trabalho contem mais cinco se§6es. A segunda se§ao realiza uma breve analise da 
literatura e resultados existentes a respeito da caracteriza§ao da evolu§ao dos pre§os 
em mercados, com o intuito de embasar a defini§ao dos modelos desenvolvidos no 
trabalho. A terceira se§ao discute prinefpios utilizados na simula§ao Monte Carlo, 
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apresentando modelos simplificados e algoritmos de simulacao utilizados nas etapas 
posteriores da pesquisa. A quarta secao envolve a caracterizacao simplificada do 
fluxo de caixa em uma operacao em mercados futuros, a partir da trajetoria de precos 
observada, encargos financeiros, tributos e taxas operacionais. A quinta secao 
apresenta os resultados da aplicacao dos procedimentos desenvolvidos nas secoes 
anteriores em estudo de caso e em modelo teorico, visando caracterizar o custo 
operacoes no mercado futuro brasileiro. Finalmente, a sexta secao apresenta 
consideracoes finais. 

2. Caracterizacao da dinamica de precos em mercados 

A dinamica dos precos em mercados vem sendo formalmente estudada pelo menos 
desde o final do seculo XIX, em larga medida com o intuito de verificar a 
possibilidade de previsao do comportamento do mercado e fundamentar a 
precificacao de derivativos. 

Contribuicao de Bachelier 

O trabalho pioneiro na caracterizacao da evolucao de precos em mercados e atribmdo 
a Louis Bachelier (1900). Em sua tese de doutorado, o matematico frances introduziu 
a nocao de que os pre§o no mercado evoluiriam de uma forma aleatoria seguindo um 
movimento Browniano, sendo, portanto, imprevisfveis. Essa pressuposi§ao teria, 
como conseqiiencia, que a esperan§a de ganhos em opera5des especulativas tenderia a 
zero. Os argumentos apresentados foram suficientes para convencer Henri Poincare, 
um dos mais famosos matematicos da epoca, que era orientador do trabalho de 
Bachelier em Sorbounne. 

Em smtese, Bachelier sugeriu, com respeito a evolu§ao dos pre§os de derivativos na 
Bolsa de Paris, que o pre§o desses derivativos num perfodo t, indicado por P t , a partir 

da informa§ao existente ate o perfodo t-l, podia ser razoavelmente representado pelo 
modelo 
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(2.1) 



onde P t e o preco no infcio do periodo t, e Ar ; e uma variavel aleatoria i.i.d. 

(independente e identicamente distribirida) ao longo dos perfodos, com E(Ar t ) = 0. O 

modelo de Bachelier assumia que as variacoes diarias Ar t tinham distribuicao normal. 

Isso era atrativo do ponto de vista matematico, mas nao eliminava a possibilidade de 
precos negativos, algo indesejavel do ponto de vista economico. 

Bachelier estudou as conseqiiencias de sua hipotese, estabelecendo bases teoricas 
para a precifica§ao de op£oes. Analisou suas ideias frente ao comportamento dos 
precos de derivativos na Bolsa de Paris, encontrando suporte empirico para o modelo 
proposto. A contribui§ao de Bachelier foi surpreendente: os metodos e resultados 
obtidos sao proximos aos existentes em desenvolvimentos recentes. Na epoca, 
contudo, o seu trabalho ficou restrito a uma pequena comunidade de matematicos - 
Kolmogorov foi influenciado por suas ideias - capaz de compreender a inovacao e 
aplicabilidade dos novos metodos desenvolvidos, argumentam Cortault et al. (2000), 
em nota biografica sobre a vida e obra de Bachelier. 

Evidencias da primeira metade do seculo XX 

Durante a primeira metade do seculo XX, desconhecendo o trabalho de Bachelier, 
autores como A. Cowles (1933), H. Working (1934) e J. Kendall (1953), utilizando 
series de pre§os de a§oes e commodities no mercado dos EUA e Inglaterra, mostraram 
evidencias fortes suportando a hipotese de imprevisibilidade dos pre§os em mercados. 
O trabalho de Cowles (1933) e particularmente interessante pelos testes que fez sobre 
qualidade das previsoes realizadas por dezenas de analistas de mercado. Apresenta 
evidencias fortes suportando a no§ao de que a sele§ao aleatoria de investimentos na 
bolsa possibilitaria chances de ganho ou prejmzo similares as sele§6es de 
investimentos realizadas por analistas profissionais. 
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No meio da decada de 50 o trabalho de Bachelier foi redescoberto por Leonard 
Savage, que o trouxe a atencao de Paul Samuelson, relata Bernstein (1992). Tanto 
Samuelson quanto seus alunos foram influenciados pelas ideias de Bachelier. Robert 
Merton, um desses alunos, viria a ganhar o premio Nobel de economia decadas 
depois, em 1997, por seu trabalho realizado na precificacao de opcoes, com modelos 
que aprimoram as ideias apresentadas originalmente por Bachelier. 

O modelo logaritmico de Osborne (1959) 

Osborne (1959) apresenta uma contribuicao importante no final da decada de 50, ao 
argumentar que seria a diferenca entre os logaritmos naturais de precos observados 
entre 2 periodos que teria distribuicao normal. O argumento de Osborne foi atrativo e 
embasado em analise da evolucao dos precos no mercado de acoes, na Bolsa de Nova 
York. Um dos aspectos interessantes do modelo apresentado por Osborne era a 
impossibilidade de precos negativos, algo que o modelo de Bachelier nao exchria. 
Esse modelo, em sua versao discreta, pode ser representado por 

ln(P f ) = ln(P f _ 1 ) + r f . (2.2) 

Nessa formulacao, r t e uma variavel aleatoria i.i.d. (independente e identicamente 

distribufda), com distribuicao Normal (0, G 2 ) . Uma conseqiiencia algebrica direta 
dessa formula§ao e 

P t = P t -ie' (2.3) 

onde e' tern distribui5ao log-normal (pelas propriedades das distribuicoes normal e 
log-normal, vide Mood et al.,1974). O termo r t , definido por \a{P t /P (-1 ) , representa a 

taxa de crescimento geometrico, em capitaliza§ao infinita. O valor de r t tende a ser, 
usualmente, muito proximo da taxa de varia§ao observada de um periodo para outro 
quando o valor da taxa e pequena (inferior a 0,10) ou seja e r ' ~ (1 + r t ) . 
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Decorre de (2.3) que o preco no periodo t tambem tera distribuicao log-normal, o que 
garante que sera estritamente positivo, algo que nao ocorrera necessariamente na 

formulacao especificada em (2.1). Em fun§ao da log-normalidade de e' , temos que 



E{e< ) = e 




(2.3.1) 



onde E( 



r t ) = \i e Var(r t ) = o 2 . Nessa situa£ao, a partir da informa§ao existente ate o 



periodo t-l, para garantir que 



E(P t ) = P t 



t-i 



(2.4) 



tal como ocorria em (2.1), e necessario que E(e') = \ ou, que a condi§ao 




(2.5) 



seja verificada. Osborne assumia que (1 = 0, o que levava a E(P t ) — P t _ { e 2 ou, 

1 , 

aproximadamente, E{P t ) = P t _ x (l +— a ) . Usava esse fato para justificar que 

tendencias historicas observadas de crescimento no pre§o de a§oes eram 
conseqiiencia do processo estocastico associado a evolu5ao dos pre§os, se (1 = 0, 
tendo pouca rela§ao com variaveis economicas. 

No modelo apresentado nos ultimos paragrafos, o desvio padrao G , e freqiientemente 
usado para caracterizar a volatilidade do comportamento dos pre§os. Nesse modelo, 
se o tempo e medido em dias, a sera caracterizado como a volatilidade diaria, que 

costuma ser expressa na base anual, por o4m , onde m e o numero de dias de 
funcionamento do mercado (245 a 260, em geral, assumido 254 nesse trabalho). Esse 
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resultado pode ser obtido facilmente, partindo do fato de que nesse modelo, o preco 
em t+m, dependera das varias taxas de variacao futuras observadas, de acordo com 



P — D „ r l+ r 2+-- + r m 

1 t+m 1 t c 



onde R m =r l + r 2 + ... + r m , representa a taxa de varia§ao no penodo m. Mas, 

V(R m )=V(r l+ r 2 + ... + r m ), 
e, pelo fato das taxas serem consideradas i.i.d e pelas propriedades da variancia, 



V(RJ = a 2 m 



Decorre dessa ultima expressao, que o desvio padrao de R m , representado por o m , e 
definido por 

o m =oVm, (2.5.1) 

onde O e a volatilidade diaria e m e o periodo considerado em dias. Essa nota§ao 
pode ser facilmente adaptada a outros perfodos que nao o diario, inclusive para 
situa§oes em que o tempo e medido de forma contmua. 



Desenvolvimentos na decada de 60 

Durante a primeira metade do seculo XX, evidencias de autocorrela§ao nas varia5oes 
de pre§os, observadas em algumas situa5oes, tendiam a ser consideradas excecao a 
regra. Working (1960) e Alexander (1961), em trabalhos independentes, mostraram 
que essas exce§6es poderiam ser explicadas por problemas metodologicos. 
Demonstram que series cujos pre§os sao medias de periodos (ex. precos mensais 
obtidos a partir da media dos pre§os diarios), podem exibir varia§6es 
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autocorrelacionadas, mesmo que as variacoes nos dados diarios nao apresentem 
autocorrelacao. 

Na decada de 60, Granger & Morgenstern (1963) e Fama (1965), utilizando um 
instrumental estatfstico mais aprimorado, realizam testes empfricos exaustivos em 
series temporais de precos, corroborando resultados obtidos anteriormente por outros 
autores. Fama (1965) argumenta que mesmo quando e observada autocorrelacao 
estatisticamente significativa, algo que mostra nao ser incomum entre variacoes de 
precos de acoes, ela teria pouco significado pratico. A autocorrelacao observada por 
Fama foi de magnitude insuficiente para justificar uma estrategia que poderia se 
beneficiar do resultado, em decorrencia dos custos de transacao existentes. Fama e 
Blume (1966) mostram outros casos em que existe alguma possibilidade de ganho 
atraves de estrategias derivadas de analises tecnicas ou grafistas, mas que esses 
ganhos tambem seriam pequenos e proximos dos custos de transacao existentes. 

Um aspecto descritivo importante analisado por Fama (1965) e o fato de que a 
distribuicao empirica das varia§6es de precos, obtida pela diferenca entre os 
logaritmos (do preco em um periodo e do preco no periodo anterior), como sugerido 
por Osborne (1959), em lugar de apresentar comportamento proximo da normalidade, 
tendia a apresentar um formato como o apresentado na Figura 2.1. Nessa figura, a 
linha cheia representa a distribuicao empirica e a linha tracejada uma normal, ambas 
com mesma esperanca e mesma variancia. A distribui§ao empirica tende a ser mais 
leptocurtica que a normal, refletindo uma curtose maior. Esse fenomeno, 
frequentemente observado na analise empirica da distribui§ao das varia§6es de pre§o, 
e caracterizado na literatura com adjetivo "cauda gorda" (fat tail em ingles). Ou seja, 
a distribui§ao tende a se parecer mais com uma t-Student (com poucos graus de 
liberdade) que propriamente com uma distribui§ao normal. Fama (1965) aponta 
algumas possfveis causas para esse fenomeno associadas a natureza do processo 
gerador das varia§oes de pre§os. Caso esse processo apresentasse variancia 
condicional nao-constante (ou mesmo altera§oes na esperan§a em certos periodos), a 
distribui§ao incondicional das varia§oes seria uma mistura de distribui§6es. O 
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resultado observado ocorreria ate mesmo se as distribui§oes condicionais fossem 
normais. 




Figura 2.1 Formato empirico da distribui§ao das taxas de varia§ao de precos (linha 

cheia) e formato da distribui§ao normal (linha tracejada), em Fama (1965) 

Samuelson (1965) e Mandelbrot (1966) interpretam as implica§oes dos resultados 
apontando para a imprevisibilidade dos precos, no contexto da teoria economica, 
sugerindo que, quando os mercados estao funcionando de uma forma apropriada, a 
informacao existente acaba sendo internalizada rapidamente pelos agentes 
economicos nos precos praticados, o que rapidamente mitigaria a possibilidade de 
ganhos especulativos. Jensen (1968) apresenta um estudo suportando essa no$ao a 
partir dos retornos observados em 115 fundos de investimento durante um periodo de 
10 anos. Mostra que as vantagens aparentes nos retornos observados em alguns 
fundos (na media essas vantagens nem mesmo existiam) acabavam se desfazendo 
quando considerados os custos de administra§ao. 

Esses resultados nao agradaram (e nao agradam) os profissionais do mercado 
defensores da analise tecnica (ou grafista) e/ou analise fundamentalista. Os grafistas 
buscam padroes na evolucao dos precos que acreditam se repetir, dando margem a 
ganhos certos em determinadas situa§6es. Os fundamentalistas, por outro lado, 
estudam exaustivamente os fundamentos dos mercados que poderiam sugerir 
estrategias vencedoras a partir da observa§ao de certos aspectos incompatfveis com os 
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precos praticados (ex. se os estoques estao baixos, entao podem existir fundamentos 
justificando o crescimento dos precos futuros, em nfveis superiores aos atuais). 

A Hipotese da Eficiencia em Mercados - Fama (1970) 

As nogoes existentes sobre imprevisibilidade foram formalizadas em Fama (1970) na 
"hipotese da eficiencia do mercado - HEM." Fama caracterizou 3 formas de 
eficiencia: fraca, semi-forte e forte. Na eficiencia fraca, toda a informa§ao passada, 
contida nos dados existentes das negocia§oes em um certo mercado, estaria refletida 
nos precos praticados, nao havendo a possibilidade de ganhos sistematicos derivados 
de estrategias fundamentadas nessas informacoes. Na eficiencia semi-forte, os precos 
praticados refletiriam, adicionalmente, toda a informa§ao dispomvel ao publico 
(incluindo a informacao dispomvel para caracterizacao da eficiencia fraca). Na 
eficiencia forte, os precos refletiriam toda a informa§ao que caracterizaria a eficiencia 
semi-forte, e, adicionalmente, refletiria tambem informacoes privadas (incluindo a 
inside information). Fama (1970) caracteriza testes para essas formas de eficiencia e 
sumariza a literatura existente ate entao. Como conclusao importante, afirma que "a 
evidencia existente em suporte da hipotese da eficiencia do mercado e extensiva e (de 
uma forma muito unica em economia) a evidencia contraria e esparsa" [traducao do 
autor]. Nos anos seguintes outros autores testaram exaustivamente as formas de 
eficiencia caracterizadas em Fama (1970), auxiliados pela crescente disponibilidade 
de novos resultados e metodologias estatfsticas, assim como recursos computacionais 
mais acessfveis. 

Modelos de precificacao de opcoes de Black- Scholes e de Merton 

Na decada de 70 ocorre um marco importante dentro dos procedimentos para 
precifica§ao de derivativos a partir de modelos que caracterizam a evolu§ao de pre§os 
de a§6es ou contratos futuros. Black e Scholes (1973), utilizando um modelo 
contmuo de varia§oes de pre§os - cujas pressuposi§6es quanta evolu5ao de pre§os 
sao comparaveis a versao apresentada por Osborne (1959) - apresentam uma solu§ao 
exata para a precifica§ao de op5oes. Robert Merton (1973) chega a resultados 
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similares, independentemente. Rapidamente o novo metodo apresentado passa a ser 
utilizado em larga escala pelo mercado, que estava carente de fundamentos mais 
solidos para precificacao de opcoes. Essa contribuicao rendeu o premio Nobel de 
economia a M. Scholes e R. Merton em 1997 (Black ja havia falecido anos antes). 

Novas metodologias dos anos 70 e 80 - Cointegra^ao, ARIMA, GARCH 

Granger e Newbold (1974) caracterizam, com um simples exemplo, uma situacao 
comum que pode condicionar uma correlacao espuria entre variaveis, que testes 
tradicionais indicariam ser estatisticamente significativa. A situacao apresentada foi a 
simples regressao linear entre observacoes temporais das variaveis Y e X, cada uma 
delas definida por um processo aleatorio como o representado em (2.1): 

Y t =Y t _,+e yt , t = l,--,n, 
X t =X t _,+e xt , t = l, — ,n, 

onde os erros observados e yt e e xt sao independentes. Numa regressao linear simples 
definida por 

Y t =b 0 +b l X t +e t , t = !,-■■, n, 

os testes estatfsticos convencionais costumam rejeitar a hipotese b x = 0 . Alem disso, 
os erros tendem a apresentar correla§ao positiva, fatos que contrariam a intui5ao. O 
uso das diferen§as entre os valores das variaveis no tempo, ou integra§ao, seria uma 
forma possfvel de tratar o problema detectado, ou seja, se a regressao for 

AF r =^AX f +e f , t = !,■•■, n, (2.6) 

onde AY, = Y t - Y t _ Y e AX t = X t - X t _ x a nulidade de b x nao sera rejeitada e os erros 
e t nao serao correlacionados. Essa observa§ao, generalizada a outras situa§6es, deu 
origem aos estudos e testes de cointegra§ao, formalizados em Engle e Grange (1987). 
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Esses estudos estabeleceram uma importante linha de pesquisa, ampliando os 
procedimentos utilizados para caracterizacao da eficiencia em mercados. 



Na area metodologica, a crescente disponibilidade de tecnicas especfficas para analise 
e caracterizacao de series temporais, durante a decada de 70, popularizaram metodos 
para identificacao e estimacao dos modelos da famflia ARIMA (AutoRegressive- 
Integration-Moving Average). Esses modelos, sistematizadas didaticamente em livros 
texto como Box e Jenkins (1976), apresentam recursos abrangentes para 
caracterizacao estatfstica de series de dados temporais. Os modelos ARIMA podem 
caracterizar, simultaneamente, num unico modelo, 3 aspectos importantes em series 
temporais: a autocorrelacao (AR), a ordem de integracao (I), e a dependencia em 
medias moveis (MA). Uma serie onde a autocorrelacao pode atingir a ordem c, 
sumarizada por AR(c), pode ser representada por 



onde e,e um erro aleatorio, i.i.d, com esperanca zero. Uma serie onde ocorrem efeitos 

associados a media movel (MA), de ordem d, sumarizada por MA(J), pode ser 
representada por 



onde e t e um erro aleatorio, i.i.d, com esperan§a zero. A situa5ao que inclui, 

simultaneamente, efeitos autoregressivos e de media movel, de ordem c e d 
respectivamente, sumarizada por ARMA(c,J), pode ser genericamente representada 
por 



Y t =Q l Y t _ l +e 2 Y t _ 2 + - + e c Y t _ c +e l 



(2.7) 



Y t = <V r-l + §2 e t-2 +■■■ + §d e t-d + 



(2.8) 



Y t = + • • • + G c F f _ c + <!>!«,_! + • • • + § d e t _ d + e, 



(2.9) 
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onde e t e urn erro aleatorio, i.i.d, com esperanca zero. O modelo pode incluir uma 
constante. 

A ordem de integracao e definida pelo numero de diferencas consideradas para a 
variavel de interesse, em que a ordem zero e a propria variavel, a ordem lea 
diferenca entre o valor emfeo valor em t-l, e assim por diante. Logo, um modelo 
ARIMA (c,i,d) indica que a autocorrelacao e considerada ate a ordem c, a integracao 
ate ordem i e o efeito de media movel ate a ordem d. Um modelo ARIMA (2,1,1) 
seria representado por: 

AY, = Q l AY t _ 1 +Q 2 AY t _ 2 + ^e,_ x +e t , (2.10) 

onde e t e um erro aleatorio, i.i.d, com esperanca zero. O modelo ARIMA basico 

apresentado, em sua versao univariada, pode ser estendido de diversas formas, com o 
objetivo de representar processos incluindo variaveis explicativas e outros aspectos. 
O diagnostico, especifica§ao e estima§ao desses modelos sao apresentados de forma 
didatica em livros texto sobre o assunto, como Hamilton (1994) e Mills (1997). 

Nos anos 80, um novo salto metodologico ocorre com o trabalho de Engle (1982). 
Esse autor - agraciado pelo premio Nobel em 2003 por essa contribui§ao - 
desenvolveu modelos estatfsticos que possibilitam caracterizar melhor certos 
fenomenos empfricos relacionados a volatilidade em series temporais de pre§os. Em 
particular, apresenta uma famflia de modelos denominada ARCH (AutoRegressive 
Conditional Heteroscedastic) para caracterizar (e estimar) o fato empiricamente 
observado de que em series de variacao de pre§os a volatilidade condicional (medida 
pela variancia) tende a ser correlacionada positivamente no tempo, apresentando 
periodos (clusters) de volatilidade mais alta e periodos (clusters) de volatilidade mais 
baixa. Bollerslev (1986) generaliza o trabalho de Engle (1982), com uma famflia de 
modelos denominada GARCH (Generalized AutoRegressive Conditional 
Heteroscedastic), a qual inclui a famflia ARCH como caso particular. 
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A representacao geral do modelo GARCH(p,q), para o caso de series univariadas, 
onde definimos por ii M a informacao dispomvel no perfodo (-1 e apresentada a 
seguir: 



Y t =E(Y t \a t _ 1 ) + ^ t 



p_ 



(2.11) 



i=\ 



onde z, e uma variavel aleatoria i.i.d.com esperanca zero e variancia unitaria. A 
distribuicao de probabilidade de z, nao precisa necessariamente ser uma normal. Os 
parametros a serem estimados sao a 0 ,---,a p e b l ,---,b g . O modelo ARCH, 
desenvolvido por Engle, seria o caso particular do GARCH, quando p=0. Nesse 
modelo, of representa a variancia condicional do erro, dada a informa§ao dispomvel 
ate o perfodo t-l, ou seja Q.,^ . Atendidas as restri§6es quanto aos valores dos 

parametros, explicitadas em Bollerslev (1986), a variancia incondicional do erro seria 
definida por 



Uma apresentacao didatica dos modelos GARCH e de outros modelos mais simples 
para tratamento da variancia condicional e apresentada por Engle (2001). Andersen et 
al. (2002, 2003) examina aprimoramentos existentes para essa classe de modelos, 
considerando generaliza§6es multivariadas, metodos alternatives de estima§ao e 
outros casos de interesse. A combina§ao dos modelos ARIMA e GARCH e uma 



0 




(2.12) 



i=i 



i=i 
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dessas possibilidades que da acesso a um amplo espectro de modelos para 
representacao do comportamento de series de precos. 

Outras estrategias complementares para estimacao de modelos flexfveis para 
caracterizacao de series temporais e previsao incluem o filtro de Kalman, cujos 
procedimentos estao sistematizados em livros texto como Harvey (1989) ou os 
modelos lineares dinamicos descritos em West e Harrison (1987,1997). 

O esforco empreendido no desenvolvimento e implementacao computacional de um 
instrumental metodologico sofisticado, especialmente a partir da decada de 70, 
possibilitou que ideias importantes sobre o comportamento dos precos nos mercados 
de acoes e mercados futuros de commodities pudessem ser verificadas. 

O massacrante escrutmio da HEM 

O poderoso arsenal de metodos desenvolvido nas decadas de 70 e 80 possibilitou uma 
melhor caracterizacao estatfstica das series temporais de variacoes de precos, 
permitindo o desenvolvimento de testes mais rigorosos para a hipotese da eficiencia 
do mercado, nos varios sentidos colocados por Fama (1970). 

Anomalias no comportamento dos precos no mercado, reveladas por diversos autores, 
estao sistematizadas em trabalhos como Fama (1991), Clarke et al. (2001), Dimson e 
Mussavian (2000) e Lima (2003). Nos mercados futuros de commodities, autores 
como Cargill e Hauser (1975), Taylor (1985) e Leuthold (1991) mostram anomalias 
que os levam a rejeitar a ideia de que modelos como os apresentados em (2.1) e (2.2) 
poderiam representar fielmente o comportamento dos precos em diversas situacoes. 

Uma linha importante de contestacao de eficiencia informacional do mercado, de 
natureza teorica, foi estabelecida por Grossman e Stiglitz (1980). Nesse trabalho, que 
sistematiza pesquisas anteriores, os autores colocam limites teoricos para a HEM a 
partir de um modelo que mostra que a validade da HEM levaria a evidentes 
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contradicts . Egelkraut et al. (2003) apresentam uma revisao da literatura e 
evidencias dentro dessa linha de pesquisa, nos mercados de milho e soja. 

Apesar das muitas anomalias reveladas e algumas dificuldades teoricas, a no§ao de 
que os precos se comportariam de uma forma imprevisfvel continua sobrevivendo 
confortavelmente como o principal modelo simplificado da realidade, usado para 
caracterizar certos aspectos fundamentais do comportamento de mercados de acoes e 
mercados futuros. Modelos de precifica§ao de op§6es como o Black-Scholes, usados 
frequentemente no mercado, assumem pressuposicoes muito simples com relacao a 
imprevisibilidade dos precos, similares as utilizadas por Osborne (1959), definidas 
em (2.3). Possivelmente, dentro de todos os pressupostos usualmente aceitos dentro 
da economia, a HEM foi o menos abalado pelos testes massacrantes aos quais foi 
submetido, concluem Clarke et al. (2001). 

Richard Roll, pesquisador reconhecido na area de finances que foi vice-presidente da 
Goldman Sachs e administrador de fundos de investimento, sugere em Roll (1994) 
uma serie de conselhos uteis para diretores financeiros de empresas. Com relacao ao 
comportamento dos precos no mercado e a possibilidade de ganhos na exploracao de 
potenciais ineficiencias, Roll argumenta que: 

"Depois de 25 anos ouvindo alegagdes sobre a ineficiencia do mercado e 10 
anos tentando explord-la eu mesmo como um gerente de investimentos, eu me 
convenci de que o conceito de eficiencia de mercado deve ser muito bem 
aprendido por diretores financeiros de empresas. A nocao de eficiencia no 
mercado e uma aproximacao de primeira ordem da realidade muito util para 
organizagao dos pensamentos sobre o funcionamento do mercado. Eu 
suspeito bastante de estudos que revelam anomalias sistemdticas ... no 
comportamento do mercado... Na ultima decada tentei explorar as 
ineficiencias mais promissoras, investindo quantias pesadas em regras de 
trading que visavam tirar proveito dessas ineficiencias... Nunca consegui 

2 Tradu^ao dos autores. 
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achar uma regra ou estrategia que funcionasse. E veja que nas fungdes que 
ocupei, freqUentemente tinha acesso as ideias mais promissoras antes mesmo 
que elas fossem publicadas ou se tornassem amplamente disponiveis, o que 
permitia testa-las em primeira mao ... O conceito de eficiencia do mercado e 
muito simples: a competigao fard com que os beneficios potenciais de uma 
nova regra de trading sejam compensados pelos custos de implementagdo ...a 
primeira pessoa que desenvolve um metodo sofisticado, uma rede neural para 
analisar o comportamento de um certo mercado, talvez consiga pagar o custo 
do desenvolvimento do metodo... mas nao a milesima. " 



Concluindo, dentro desse trabalho, sao apresentados procedimentos que possibilitam 
a simulacao de trajetorias de precos tanto com modelos simples como o descrito em 
(2.2) e (2.3), como por modelos mais complexos e mais aderentes a realidade, 
derivados da implementacao de modelos ARIMA e GARCH. Uma preocupacao 
importante e a verificacao do impacto do grau de refinamento utilizado na solucao 
dos problemas tratados na pesquisa. 



3. Simulacao de trajetorias de precos e funcoes destas pelo Metodo Monte Carlo 

O uso da simula§ao em modelagem e uma tecnica fundamental para a solu§ao 
problemas quando e diffcil, muito trabalhoso, ou mesmo impossivel, o uso de 
metodos analfticos exatos. A modelagem atraves de metodos analfticos exige modelos 
parcimoniosos que podem ser pouco aderentes a realidade em muitos casos. Atraves 
de tecnicas de simula§ao, e possfvel o desenvolvimento de modelos complexos, mais 
aderentes a realidade. Herbert Simon (1996) caracteriza a simula§ao como uma 
tecnica essencial para o entendimento do funcionamento de sistemas complexos, 
especialmente quando o comportamento desse sistema nao pode ser facilmente 
inferido a partir das propriedades das partes que dele fazem parte. Em muitas 
situacoes e diffcil estabelecer as implica§6es condicionadas por certas premissas. 
Nesses casos, a simulacao pode ser util para caracterizar as implica§6es derivadas das 
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premissas consideradas, gerando com isso conhecimentos previamente inexistentes. 
Dentro deste trabalho, o problema central e justamente o entendimento das 
conseqiiencias de pressuposicoes alternativas sobre o comportamento das trajetorias 
de precos em mercados futuros sobre fluxos de caixa e custos operacionais, algo 
pouco formalizado. 

No passado, a utiliza§ao da simula§ao era dificultada pela disponibilidade limitada de 
recursos computacionais. Especialmente quando eram necessarias muitas simulacoes 
para obtencao dos resultados de interesse dentro de limites preestabelecidos de 
precisao, num intervalo de tempo aceitavel. A disponibilidade cada vez maior de 
recursos computacionais, ao longo das ultimas duas decadas, vem favorecendo a 
utilizacao de solu§6es fundamentadas em simulacao. 

A simula§ao estocastica ou realizada a partir do Metodo Monte Carlo, na sua forma 
mais basica, parte de um procedimento (um gerador) que possibilita o sorteio de 
numeros aleatorios entre 0 e 1, para sortear numeros com qualquer distribuicao de 
probabilidade. 

A base teorica que garante esse resultado e o teorema da transformacao integral de 
probabilidades, detalhado em Mood et al. (1974). Esse teorema garante que, a partir 
de uma seqiiencia de numeros aleatorios entre 0 e 1, com distribuicao uniforme e 
independente, e da funcao inversa da fun§ao de distribuicao probabilfstica que 
caracteriza o fenomeno que desejamos simular, podemos gerar uma seqiiencia de 
numeros aleatorios com essa distribui§ao desejada. Ocorre, contudo, que muitas vezes 
e dificil inverter, analiticamente, a fun§ao de distribui§ao de interesse. Nesses casos 
pode-se recorrer a uma aproxima§ao dessa funcao, ou mesmo utilizar outro metodo 
de simulacao. Bratley et al.(1983), Law e Kelton, (1991) e Robert e Casella (1999) 
descrevem diversas tecnicas para a simulacao de numeros aleatorios tais como a 
aceita§ao-rejei§ao ou metodos especfficos envolvendo transforma§ao, assim como 
procedimentos que visam otimizar a simula§ao. 
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A geracao de numeros aleatorios entre 0 e 1 com distribuicao uniforme, 
independentes, algo aparentemente simples, nao e tarefa trivial. Os geradores 
existentes, implementados em software, produzem seqiiencias que se aproximam 
bastante do ideal, mas tern limitacoes, sendo por isso chamados de geradores de 
numeros pseudo-aleatorios. 

A simulacao Monte Carlo (e suas extensoes) e o metodo mais geral disponfvel para 
caracterizacao estatfstica de funcoes de variaveis aleatorias. No presente trabalho, a 
tecnica e usada exatamente nesse contexto, para conhecer estatfsticas de funcoes nao 
triviais de uma trajetoria futura de precos de um dado contrato no mercado futuro. 
Para melhor caracterizar o problema, assumimos que estamos no perfodo 0 e usamos 

P.AA-^m (3-D 

para representar as variaveis aleatorias que definem os precos de ajuste ao final de 
cada perfodo desde o perfodo 1 ate perfodo m, quando termina a operacao. Suponha 
que o interesse presente e sobre a distribui§ao de probabilidade de uma variavel 
aleatoria Y, definida como sendo uma fun§ao g(.) das variaveis aleatorias que definem 
a seqiiencia de pre§os futuros ou seja 

Y = g(P l ,P 2 ,P 3 ,-;PJ. (3.2) 

Na apresenta5ao a seguir, a fun§ao g(.) pode ser univariada (mxl) ou, em geral, 
multivariada (mXv) , onde v e a dimensao de Y para cada trajetoria de pre§os. A 
fun§ao g(.) poderia ser uma fun§ao de interesse, por exemplo, max(P 1 ,P 2 ,P 3 ,---,P m ) , 

ou seja, o maximo valor observado nos m perfodos. A fun§ao pode depender de 
outros parametros, como na fun§ao max(K-P m ,0), o resultado de uma opcao put 

europeia onde K e o pre50 de exercfcio. No caso multivariado mx2 a funcao g(.) 
poderia ser, por exemplo, 
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g(p 1 ,p 2 ,p 3 ,-,PJ= 



max(P l ,P 2 ,P 3 ,--;PJ 
max(£-P m ,0) 



(3.2.1) 



Na proxima secao caracterizaremos funcoes relativamente nao triviais, associadas ao 
fluxo de caixa de operacoes em mercados futuros. 

Se a distribuicao de probabilidade multivariada de P l ,P 2 ,P 3 ,---,P m fosse 

perfeitamente conhecida e fosse disponfvel um metodo para sortear (usando tecnicas 
de simulacao) seqiiencias a partir dessa distribuicao de probabilidade, poderiamos 
obter o valor de Y resultante de cada uma dessas seqiiencias. Supondo que obtivemos 
n seqiiencias por 



onde o primeiro subescrito representa o numero da seqiiencia simulada, e possfvel 
responder qualquer questao sobre a distribuicao de Y, lan5ando-se mao da Lei dos 
Grandes Numeros, um dos principals resultados da estatfstica. No caso de g(.) ser 
uma fun§ao univariada, a esperan§a e a variancia, assim como todos os momentos da 
distribui§ao de probabilidade de Y, , podem em geral ser facilmente estimados pelas 
contrapartes amostrais: media, variancia amostral e momentos amostrais, em 
decorrencia da Lei dos Grandes Numeros: 



ou seja, o valor do momento absoluto amostral de ordem k converge (em 
probabilidade) para o momento absoluto teorico de ordem k, para a variavel aleatoria 
Y. Esse resultado somente depende de condi§6es de regularidade como a existencia 
dos momentos teoricos (Mood et al., 1974). Por exemplo, a esperan§a matematica que 






(3.3) 



plim^^ = E(Y k ) 



(3.4) 
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corresponde ao momento absoluto teorico de ordem k=l, pode ser estimada pelo 
momento absoluto amostral de ordem k=l. A precisao dos momentos amostrais, 
encontrados a partir de n simula§oes, pode estimada atraves do intervalo de confian§a 
apropriado. Para a esperanca matematica, por exemplo, o intervalo de confian§a, a um 
nfvel de confian£a p, e estimado por 

Y±t (3.4) 

n — p,n-\ / v / 

■sjn 

onde Y n 6 a. media, s n 6 o desvio padrao amostral, obtidos a partir de n simulacoes, e 
t p ,n-\ & 0 valor critico que define, para uma variavel aleatoria T com distribui5ao t- 
Student com n-1 graus de liberdade, 

Pr ( r< v „_ 1 ) = f • (3.5) 
Como n e usualmente um numero superior a 1000, t p n _ x pode ser bem aproximado 
por z p , o valor critico que define, para uma variavel Z com distribui5ao normal 
padronizada, 

Vx{Z<z p ) = ^. (3.6) 

Os valores de z p para mveis criticos usuais, considerando nfveis de significancia 

correspondentes a p=0,95 e p=0,99, sao, respectivamente, 1,960 e 2,576 
(aproximados para 3 decimals). 

Frequentemente, a simula§ao e usada para estimar 

Pr(a < Y < b) (3.7) 

ou seja, a probabilidade de observarmos realiza5des da variavel de interesse Y dentro 
do intervalo [a,b]- Nesse caso, basta que seja definida uma variavel aleatoria X, que 
assume 1 quando a<Y <b e valor 0 em outros casos. 
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Do ponto de vista teorico, como essa variavel X, por definicao, tera distribuicao de 
probabilidade Bernoulli, sera verdade queE(X) = Pr(a < Y < b) , dado que os eventos 
a <Y < b , e X =1 sao equivalentes do ponto de vista probabilfstico. Esse fato leva a 
probabilidade de X = 1 a ser identica a Pr(a <Y <b) . A estimativa de 

Pr(a < Y < b) pode, entao, ser obtida pela media de X v X 2 , • ■ • , X m , que nada mais e 
do que a freqiiencia de observa§6es onde X = 1 . 

Os valores de X 1 ,X 2 ,---,X n podem ser encontrados a partir dos valores de 
Y 1 ,Y 2 ,---,Y n , simulados a partir dos procedimentos explicitados em (3.5). Nesse caso, 

a precisao pode ser estimada a partir do intervalo de confianca derivado anteriormente 
em (3.6), ou mesmo a partir de metodos estatfsticos Bayesianos. Nesse ultimo caso, a 
distribui§ao exata de probabilidade de Pr(a < Y < b) , derivada a partir da n 
simula§6es, pode ser representada uma distribui§ao Beta, com parametros definidos a 
partir do numero de observa§oes em que X = 1 . Esse resultado depende da 
pressuposi§ao de que o conhecimento a priori sobre Pr(a < Y < b) e representado por 
uma distribui§ao uma distribui5ao Uniforme (0,1). 

O proprio formato da distribuicao de probabilidade de Y pode ser estimado 
diretamente dos valores simulados, a partir de procedimentos descritos em Silverman 
(1989) e Scott (1992). 

Os resultados apresentados para o caso de g{.) univariada podem ser generalizados 
para o caso multivariado mXv . Nessa situa§ao, cada trajetoria de pre§os resultara 
num vetor Y de dimensao v > 1 . Tambem nesse caso, a Lei dos Grandes Numeros 
possibilita a estimativa do grau de associa§ao estatfstica dos elementos desse vetor 
atraves das estatfsticas amostrais. Em particular, se v = 2 , seria possfvel utilizar a 
correla§ao amostral entre as 2 variaveis que caracterizam o vetor Y , a partir das n 
simulacoes realizadas. Essa correlacao amostral, pela Lei dos Grandes Numeros, e 
atendidas condi§6es de regularidade, converge em probabilidade para a correla§ao 
teorica. 
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Exemplos e Algoritmos 

A seguir sao apresentados 3 exemplos que ilustram procedimentos descritos nos 
paragrafos anteriores, caracterizados por algoritmos simplificados, formalizados em 
linguagem estruturada. Sao desenvolvidos em grau crescente de complexidade, com o 
objetivo de facilitar o entendimento de um algoritmo geral para simulacao de series 
de precos definidas a partir de um processo ARIMA-GARCH, e dos 
desenvolvimentos das proximas secoes. Ilustram tambem a adicao de mais 
complexidade aos modelos visando uma melhor aderencia a realidade e as possfveis 
conseqiiencias em termos de custos e beneffcios. Para cada exemplo, e apresentado 
um algoritmo para obtencao da solu§ao desejada. As situa§oes tratadas incluem 
diferentes modelos estatfsticos usados para caracterizacao das series temporais de 
precos, assim como problemas tecnicos associados a sua implementacao. Esses 
algoritmos foram, posteriormente, implementados em programas desenvolvidos na 
linguagem associada ao software estatfstico R (versao 1.7.1), para efeito de obtencao 
dos resultados desejados. 



Exemplo A - Um investidor, ao final do dia 0, quando o preco de fechamento 
observado para o contrato futuro de certa commodity foi $ 15, deseja realizar 
inferencias sobre o maximo valor que pode atingir o preco de fechamento desse 
contrato durante os proximos 50 dias de operacao desse mercado. Como esta 
fazendo uma operacao de hedge, esse conhecimento pode ser util para 
dimensionamento dos recursos necessarios para fazer frente a potenciais ajustes 
diarios. Estudos empfricos sobre esse contrato indicam que o comportamento dos 
pre§os pode ser representado a partir do modelo sugerido originalmente por 
Osborne (1959), formalizado anteriormente pelas equa§oes (2.2) a (2.5). Os 
estudos empfricos realizados estimaram o valor da volatilidade diaria, O , como 
sendo igual a 0,016, correspondente a 0,25, aproximadamente, em base anual 
(considerando 254 dias de opera§ao no ano). O investidor assume que o pre§o no 
dia 0 e uma estimativa nao viesada do pre§o ao final dos 50 dias. 
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Para resolver esse problema e necessario que se realize, inicialmente, a simulacao de 
n trajetorias de 50 precos a partir da distribui§ao multivariada de probabilidade de 



P\, P 2 > °3' " ' ^50 (3-8) 

que depende, nesse caso, das pressuposicoes descritas em (3.1) - (3.3) e do fato 
conhecido, P Q = 15. Para definicao da trajetoria i=l, das n trajetorias que desejamos 
simular, podemos simular 

r n> r i2> r i3> — ,r iiS0 , i = l, (3.8.1) 

onde r~6 um valor simulado de uma distribuicao Normal(-^-o 2 ,o 2 ), onde 

G = 0,016 e a volatilidade diaria especificada no problema. Considera-se p. = -^-a 2 , 

para incorporar a nocao de que o pre§o no periodo 0 e uma estimativa nao 
tendenciosa do pre§o no periodo m, conforme a discussao apresentada na secao 2 do 
trabalho. Outras pressuposi§oes poderiam ser utilizadas. 

A partir dos valores simulados, seria definida a primeira trajetoria de pre§os pelo 
procedimento recursivo definido por 

P Ut =P Ut _/\ f = l e ? = l,---,50, (3.9) 

onde P i<0 = 15, o que definiria a trajetoria 1 por 

P P P ■■■ P 

1 W 1 12' ^ 13' '' l l,50- 

Para essa primeira trajetoria, o valor de 5^ , resultado derivado da fun§ao da trajetoria 
de interesse, genericamente especificada por g(.) em (3.2) e (3.3), e o maximo dos 
valores observados em uma dada trajetoria. Para i=\, temos 
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^=max(^,^. 2 ,^3,---,^5o). 



(3.10) 



Para simulacao das outras trajetorias, de 2 a n, o processo realizado em (3.8) a (3.10) 
deve ser repetido com i=2, ...,n, o que possibilitaria a obtencao da seqiiencia de n 
valores maximos ou seja 

Y x ,Y 2 ,Y„-,Y n (3.11) 
Os procedimentos apresentados nos ultimos paragrafos podem ser sumarizados pelo 
algoritmo 3A descrito na Figura 3.1, onde g(.) e a funcao max(.), m=50 e P 0 = 15. 

No algoritmo, os indices correspondentes as simulacao foram omitidos em algumas 
variaveis, em funcao de nao serem de fato necessarios para o processo computacional. 
Em cada simula§ao (e as vezes em cada periodo) as mesmas variaveis sao redefinidas 
sempre que possfvel, sem qualquer prejufzo. O programa que implementou o 
algoritmo foi utilizado tambem para estimar a esperan§a, o desvio padrao, a 
probabilidade do maximo preco ser igual ou superior a 20 e a propria distribui§ao de 
probabilidade de Y, a partir dos valores simulados e procedimentos descritos em 
paragrafos anteriores. 

As estimativas foram obtidas a partir de experimentos que consideraram 1.000, 
10.000 e 100.000 simulacoes, representadas nos algoritmo por n. A Tabela 3.1 
apresenta o resultado dessas simula§6es. A Figura 3.2 apresenta uma estimativa da 
propria distribui§ao de probabilidade, a partir de procedimentos descritos em 
Silverman (1989) e Scott (1992), considerando um kernel gaussiano, utilizando o 
experimento com 100.000 simula§6es. 
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Algoritmo 3A 

inicio 

especifique: P 0 , O, \i,m, n, fungaog(.) 

para i <— 1, • • • , n , faga { 

para t <— 1, • • • , m faga { 

simule r~ Normal(\i,<5 2 ) 

limite r t a maxima oscilagao permitida 

} 

Y t ^g(P lt P 2t P 3t - t P H ) 

} 

estime momentos amostrais a partir de Y l , Y 2 , Y 3 , • • • , Y t 

estime probabilidades de eventos de interesse 
estime a distribuigao de probabilidade de Y 

fim 



m 



Figura 3.1 Algoritmo utilizado para simula§ao de funcoes de trajetorias de pre50s, 
a partir das taxas de variacao (Exemplo A). 
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Parametro 


Estimativas obtidas em cada 
Experimento (numero de simulacoes) 


1.000 


10.000 


100.000 


E(Y) 

Intervalo confianca a 99% 


16,173 
16,08 a 16,26 


16,243 
16,21 a 16,27 


16,242 
16,23 a 16,25 


Pr(y > 20) 
Intervalo confianca a 99% 


0,0060 
-0,0003 a 0,0012 


0,0076 
0,0054 a 0,0098 


0,0076 
0,0069 a 0,0083 



Tabela 3.1 Estimativas da E(F) e Pr(Y > 20) pela media amostral e seus intervalos 
de confian§a a 99% de probabilidade (Exemplo A, Algoritmo 3A) 




Figura 3.2 Estimativas da distribui§ao de probabilidade de Y (Exemplo A, Algoritmo A) 
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Exemplo B - Mesma situacao do Exemplo A, com a diferenca que nesse caso 
observa-se empiricamente que as taxas de variacao dos precos tern uma 
autocorrelacao de ordem 1 . Especificamente, 

rt=Q r t-i+$t> (3-12) 

£ ~ Normal^, 0 2 ). (3.13) 
Esse modelo e identico ao utilizado no Exemplo A, quando 9 = 0. Deseja-se 

estimar a esperanca do maximo valor observado, para 0 = 0, 0 = 0,10, 0 = 0,20 
e 9 = 0,30 e ties mveis de volatilidade diaria, 0,0125, 0,0160 e 0,019, que 
correspondem, aproximadamente, as volatilidades anualizadas 20%, 35% e 30%. 
Os experimentos deverao ser realizados com 10.000 simulacoes. Tem-se tambem 

como informacao conhecida o valor de r 0 , 0,002, obtido por r 0 = ln(P 0 / P_ { ) , a 
partir de informa5des dispomveis no periodo 0. Considere p.=-^-o 2 , como no 
exemplo anterior. 

A autocorrela§ao existente nesse modelo traz algumas complica§6es que devem ser 
examinadas. Em geral, na simula§ao, em decorrencia de (3.12) e condicional a 
informa§ao existente no tempo t-l, temos 

E(r t ) = Qr t _ l + E(^) (3.14) 
V(r t ) = a 2 (3.15) 

Como 

E(P t ) = P t _ x E(e r >) 

usando o resultado teorico de (2.3.1), e pressuposi§oes utilizadas no Exemplo A, chega-se 
a 



E(P t ) = P t _/^ (3.16) 
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ou seja, somente se 9 = 0 fica garantido E(P t ) = P t _ x independentemente dos valores 
observados para r t . 

O algoritmo proposto para solucao do problema especificado no Exemplo B, para 9 e O 
(volatilidade diaria) definidos, e apresentado na Figura 3.3. Os resultados dos 
experimentos requeridos para resposta as questoes de interesse, obtidos a partir de 10.000 
simulacoes, sao apresentados nas Tabelas 3.2 e 3.3. Esses resultados mostram o impacto 
de introduzir a autocorrelacao nas a trajetarias de precos e na estimativa da esperanca do 
maximo valor observado, assim como na probabilidade de termos em algum momento 
precos futuros iguais ou superiores a $ 20. O efeito de um 0 mais alto tende a ser 
amplificado em cenarios de volatilidades mais elevadas. Tanto 0 quanta (J contribuem 
positivamente para o crescimento da esperanca e da probabilidade de valores mais 
extremos, um resultado que tambem pode ser observado teoricamente, a custa de 
desenvolvimentos analfticos. 
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Algoritmo 3B - AR(1) 

inicio 

especifique: P 0 , r 0 , 0, [i,o , m, n, fungaog(.) 
para i <— 1, • • • , n , faga { 

para t <— 1, • • • , m faga { 

simule £ ~ Normal ( (i , o 2 ) 

limite r t a maxima oscilagao permitida 

i 

Y t ^g(P lt P 2t P 3t - t P m ) 

} 

estime momentos amostrais a partir de Y 1 , Y 2 , Y 3 , • • • , Y n 

estime probabilidades de eventos de interesse 
estime a distribuigao de probabilidade de Y 



Figura 3.3 Algoritmo utilizado para simula§ao de fun§6es de trajetorias de pre§os, 
com autocorrela5ao das taxas de ordem 1 - AR(1), (Exemplo B). 
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A 

u 


Volatilidade Anualizada 


0,20 


0,25 


0,30 


0 


15,97 


16,25 


16,48 


0,10 


16,06 


16,38 


16,61 


0,20 


16,18 


16,51 


16,82 


0,30 


16,33 


16,74 


17,05 



Tabela 3.2 Estimativas da E(Y) pela media amostral (Exemplo B, Algoritmo 3B) 



e 


Volatilidade Anualizada 


0,20 


0,25 


0,30 


0 


0,001 


0,007 


0,021 


0,10 


0,002 


0,018 


0,038 


0,20 


0,007 


0,028 


0,062 


0,30 


0,012 


0,048 


0,089 



Tabela 3.3 Estimativas Pr(Y > 20) pela media amostral (Exemplo B, Algoritmo 3B) 



Exemplo C - Mesma situacao do Exemplo A, com a diferenca que nesse caso 
observa-se empiricamente que as taxas de variacao dos precos tern uma 
autocorrelacao de ordem 1 e que a variancia condicional do erro nao e 
homogenea, comportando-se como num modelo GARCH(1,1). Essa especificacao 
e comum como representacao adequada de muitas series temporais de variacao de 
precos. Especificamente, 

r t =0r M ,+$ ( (3.17) 




(3.18) 



6> t =a 0 +atf_ l +h&_ l (3.19) 

Z, ~Normal(0,l) (3.20) 

Esse modelo e identico ao utilizado no Exemplo A, quando 
0 = 0, a { =0, e Z?j=0. Nas expressoes (3.18) e (3.19), <5 2 representa a 

variancia incondicional e G 2 representa a variancia condicional a informacao 
existente no periodo t-\. Nesse caso deseja-se estimar a esperan5a do maximo 
valor observado dos pre5os, para 9 = 0, 0 = 0,10, 0 = 0,20 e 0 = 0,30, e para 
tres nfveis de volatilidade diaria, incondicionais, de 0,0125, 0,0160 e 0,019, que 
correspondem, aproximadamente, as volatilidades anuais, incondicionais, 0,20, 
0,25 e 0,30. Os experimentos deverao ser realizados com 10.000 simulacoes. 
Tem-se tambem como informa§ao conhecida o valor de r 0 , 0,002, obtido por 
r 0 =ln(P 0 /P_ l ) , a partir de informacoes disponfveis no periodo 0. Os valores dos 

parametros que definem a variancia condicional sao a x = 0,3 e b x = 0,2 . Os 
valores de fl 0 considerados sao, respectivamente, 0,00007825, 0,000128 e 

0,0001805, os quais levam o desvio padrao incondicional aos valores 0,0125, 
0,0160 e 0,019, em razao de (2.12). A variancia condicional, no periodo 0, e 
assumida identica a variancia incondicional (em outras situacoes poderia ser 
representada pela variancia incondicional estimada para o periodo 0). 
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O algoritmo proposto para solucionar o problema e apresentado na Figura 3.4. Esse 
algoritmo foi implementado na linguagem R, segundo procedimentos ja descritos. Os 
resultados obtidos sao similares aos observados nas Tabela 3.2 e 3.3. Para o caso 
extremo tratado no Exemplo B, onde 0 = 0,30 e a volatilidade anualizada e 0,30, um 
experimento com 10.000 simulacoes encontrou Y =17,06, com intervalo de 
confianca para a 99%, definido por [17,01; 17,11]. O valor encontrado para o 
exemplo B, nas mesmas condicoes com relacao a volatilidade e autocorrelacao, vide 

tabela 3.2, foi Y = 17,05 . As estimativas de Pr(Y > 20) , obtidas para os casos B e C, 
foram tambem muito proximas. Esse exemplo mostra que, no caso da funcao g(.) 
considerada, a modelagem da dinamica dos precos de uma forma mais complexa nao 
foi capaz de modificar a magnitude das estimativas desejadas e introduziu custos e 
dificuldades desnecessarias. O experimento que utilizou o algoritmo C demorou cerca 
de 3 vezes mais tempo para ser executado que o realizado com o algoritmo B. 
Adicionalmente, a maior complexidade do algoritmo C aumenta as chances de erros 
nas v arias etapas requeridas para sua implementacao. 

Os exemplos anteriores introduziram procedimentos de simulacao que inclufam, em 
ordem crescente de complexidade, aspectos empfricos freqiientemente observados em 
series temporais de variacoes de precos. Para completar essa secao, apresentamos um 
algoritmo mais geral, na Figura 3.5, voltado a simulacao de series de precos cujas 
taxas de variacao podem ser representadas por um processo ARIMA(c,0,d)- 
GARCH(p,q). Na situacao mais elementar, quando os parametros associados aos 
processos ARIMA e GARCH sao identicos a zero, o algoritmo e funcionalmente 
equivalente ao algoritmo A. O algoritmo 3D depende do conhecimento inicial de uma 
serie de parametros, usualmente estimados a partir de procedimentos apropriados. Em 
particular, depende do conhecimento da ordem e parametros dos modelos ARIMA e 
GARCH, do pre§o no perfodo 0, das variancias condicionais de periodos anteriores a 
0, da esperan§a de \ t , representada por [i , que pode ser especificada por (2.5) ou 

outra pressuposi§ao, das taxas r t e valores de \ t observados antes e ate perfodo 0, 

para atendimento a inicializa§ao do modelo. A defini§ao dos valores de \ t 
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observados antes do penodo 0, a partir de (X e o (a variancia incondicional), pode 
ser realizada a partir do procedimento pratico proposto, por exemplo, em Laurent e 
Peters (2002), o qual e considerado no algoritmo 3D, o que nao exclui outras 
possibilidades, inclusive a utilizacao de resfduos de modelos estatfsticos 
desenvolvidos para estimacao dos parametros. 



Algoritmo 3C - AR(1)-GARCH(1,1) 

inicio 

especifique: P 0 , r 0 , 9, n,c , c 0 , a 0 , a x , b x , n, m, fungaog(.) 
para i <— 1, • • • , n faga { 

simule % 0 ~ Normal ( \i,G 2 ) 

para t <— 1, • • • , m faga { 

simule z ~ Normal(0,l ) 



V a o+«i°f-i+^-i 



t 1 2 

limite r t a maxima oscilagao permitida 

} 

Y t ^g(P lt P 2 ,P 3 ,:;P m ) 

} 

estime momentos amostrais a partir de Y l ,Y 2 ,Y 3 ,---,Y n 

estime probabilidades de eventos de interesse 
estime a distribuigao de probabilidade de Y 

fim 



Figura 3.4 Algoritmo utilizado para simula§ao de fun§6es de trajetorias de pre50s, 
com autocorrela5ao das taxas de ordem 1, AR(1) e variancia condicional 
nao homogenea(Situa§ao C) representada por um processo 
GARCH(U) 
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Algoritmo 3D - ARIMA(c,0,d)-GARCH(p,q) 

inicio 

especifique: P 0 , a , \i , a Q , n, m, p, q, c, d, fungaog(.) 

a l ,a 2 ,---,a p ,b l ,b 2 ,---,b q 

r o ' r -i ' " ' ' r l-c 
O0' a -l'""'°l-p 

para i <— 1, • • • , n faga { 

simule £o > 4-i > • ' ■ki-m^d^) ~ Normal( y.,0 2 ) (ou inicializado) 
para t <— 1, • • • , m faga { 

simule z ~ Distribuigao( 0,1) (Normal, t-Student, etc.) 

(veja2.11) 

//m/fe r ( a maxima oscilagao permitida 

} 

estime momentos amostrais a partir de Y l , Y 2 , Y 3 , • • • , Y n 

estime probabilidades de eventos de interesse 
estime a distribuigao de probabilidade de Y 



Figura 3.5 Algoritmo utilizado para simulacao de funcoes de trajetorias de precos, 
com autocorrelacao das taxas de ordem 1, AR(1) e variancia condicional 
nao homogenea(Situacao C) representada por um processo 
GARCH(U) 
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4. Caracterizacao do fluxo de caixa e custos operacionais em mercados futuros 

Esta secao examina o problema da caracterizacao do fluxo de caixa de um agente 
economico em operacoes em mercados futuros. Um dos objetivos dessa analise e o 
interesse existente pela modelagem e caracterizacao de custos operacionais, decorrentes 
de margem de garantia, ajustes diarios, taxas operacionais e tributos. 

A analise assume que o tempo e medido em dias de operacao do mercado e que o agente 
economico inicia a operacao no periodo 0, adquirindo um contrato futuro de venda que 
ira liquidar no periodo m. O resultado dessa operacao, no periodo m, pode ser 
representado por 

R = A-CT(P Q f ,P/,-;Pj l ,K) (4.1) 

onde A e o valor do ajuste acumulado, definido por Pq - P„ , CT(Pq , P/ , • • •, P^ , 
sao custos de opera§ao, dependentes da trajetoria dos precos observada no mercado 

futuro, e de K, um vetor de parametros que caracterizam taxas de juros, tributos, taxas 
operacionais etc., assumidos como conhecidos no infcio da opera§ao. O pre§o utilizado 
reflete o valor total do contrato, e nao o pre§o unitario da commodity considerada nesse 
contrato. Essa nota§ao visa facilitar a analise, dado que alguns custos sao especificados 
por contrato. Como as estimativas de custo serao expressas no final da analise em termos 
percentuais, e possfvel traduzir facilmente os custos e seus varios componentes, em bases 
comparaveis. 

Em geral, alguns desses parametros podem ser incertos (ex. taxas de juros), alterando-se 
ao longo da opera§ao, mas serao considerados conhecidos e fixos para efeito deste 
trabalho. A incerteza com rela§ao ao valor desses custos vem da incerteza associada as 
possfveis trajetorias que o pre§o futuro pode seguir entre os perfodos 0 e m. Cada 
trajetoria define um unico custo e por isso podemos dizer que o custo e uma funcao da 
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trajetaria e usar, para caracterizar esse custo, os procedimentos desenvolvidos na secao 3, 
a partir da simulacao de trajetarias de precos segundo um modelo apropriado. 

Dentro do modelo descrito a seguir, alguns aspectos de operacoes reais sao abstrafdos, 
sem grande prejufzo, com o objetivo de facilitar a exposicao. Os aspectos praticos 
abstrafdos podem ser implementados sem grandes dificuldades em versoes mais 
aprimoradas do modelo. 

Serao considerados 5 tipos de custo: (a) taxas operacionais (taxa de registro, taxa 
operacional basica (TOB) e emolumentos); (b) desembolsos associados ao IRPF; (c) 
desembolsos associados a CPMF; (d) encargos financeiros referentes a margem de 
garantia; e, (d) encargos financeiros associados aos ajustes diarios. A situacao modelada e 
a de um agente do tipo pessoa ffsica. Para pessoas jurfdicas, tributadas com base no lucro 
real ou no lucro presumido, instituicoes isentas de certos tributos etc., os custos 
operacionais observados podem ser diferentes dos apresentados neste trabalho, sendo 
assunto de pesquisas futuras. 

Taxas operacionais e tributos 

As taxas operacionais incluem 3 componentes: (a) taxa de registro por contrato, ex. R$ 
0,20 para contratos de commodities agricolas negociados em R$ e dolar; (b) TOB (taxa 
operacional basica), correspondente a 0,3% do valor do contrato, para commodities 
agricolas e dolar, por exemplo; (c) emolumentos, calculado como um percentual de 
6,32% da TOB (o valor e negociavel com a corretora em funcao do volume). A taxa de 
liquidacao por entrega sera desconsiderada, assumindo-se que o encerramento da 
operacao sera realizado por reversao da posicao. As 3 taxas sao cobradas tanto no initio 
quanta no final da operacao. A TOB e calculada sobre o preco de ajuste do dia anterior 
do contrato correspondente ao segundo vencimento. Para simplificar a aplicacao dessas 
taxas, o trabalho assume que a TOB sera calculada sobre o preco do contrato no periodo 

3 Os dados utilizados sobre tributos e taxas operacionais, descritos nesta secao foram obtidos dentro da 
pesquisa realizada junto a BM&F e agentes do mercado, conduzida dentro da dissertacao de mestrado de E. 
Andrade. 
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0. Os valores apresentados referem-se a operacoes normais, pois operacoes de day trade 
consideram taxas diferentes das apresentadas (veja no site da BM&F). 

O imposto de renda (IRPF) nas operacoes em mercados futuros e regulamentado pela IN 
da SRF 25/01. As regras para calculo dos valores devidos sao descritas a seguir. O IRPF, 
na alfquota de 20%, e cobrado dos saldos positivos resultantes da soma algebrica dos 
ajustes diarios, apurados ao final de cada mes do calendario. Se esse saldo for negativo 
num dado mes, esse valor sera transferido para o calculo do saldo nos meses 
subseqiientes. O IRPF deve ser pago 15 dias uteis apos o final do mes da apuracao do 
tributo. Considera-se, contudo, neste trabalho, para simplificar o modelo, e sem 
alteracoes sensfveis nos resultados, que o imposto e pago ao final do mes dentro qual ele 
e devido, sendo o valor descontado pela taxa de juros correspondente ao custo de 
oportunidade, o que considera uma antecipacao de pagamento por 15 dias. 

No caso da CPMF, considera-se o valor da alfquota como sendo 0,38%, sendo cobrada 
(de acordo com a lei 9311/96) sobre: (a) os pagamentos das taxas operacionais; (b) 
deposito de margens de garantia em dinheiro; (c) o saldo lfquido de ajustes diarios de 
toda a operacao, se esse saldo for negativo. Nao sao consideradas outras tributacoes por 
CPMF que poderiam ocorrer em outras movimentacoes financeiras utilizadas para a 
gestao da operacao, o que pode subestimar os desembolsos para pagamento desse tributo. 

Custos financeiros associados a margem de garantia e aos ajustes diarios 

Na BM&F, a margem e definida como um valor fixo por contrato, que pode ser 
aumentado por chamadas de margem, durante a operacao, em condicoes especificas do 
mercado. Para o contrato de milho, por exemplo, que corresponde a 450 sacas de 60 kg, o 
valor da margem era R$ 560,00, em 30 de outubro de 2003, o que representava 5,6 % a 
6,5 % do valor dos contratos nessa data. Para operacoes fora do hedge, a margem era R$ 
700, o que representava entre 7 e 8 % do valor do contrato. A margem pode ser 
depositada em certificados de deposito bancario (CDBs) e por outros ativos. Para o 
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contrato de dolar comercial, essa margem variava, na mesma data, entre R$ 13.731,00 
(NOV 03) e R$ 34.898,91 (OUT 94), para um contrato de US$ 50.000,00. 

Para formalizacao dos custos financeiros, assume-se que o agente tern pela frente duas 
taxas de juros relevantes: (a) uma taxa J a que remunera suas aplicacoes, refletindo o 
custo de oportunidade do seu capital, pressuposta identica a taxa media de remuneracao 
de CDBs; e (b) uma taxa J c para captacao de recursos. A ideia geral do modelo 

apresentado e que o agente recebe juros J a pelos saldos positivos de uma conta que 
reflete a operacao e paga juros J c pelos saldos negativos. O recebimento de juros sobre 
saldos positivos e facilitado por mecanismos como o FIF (Fundo de Investimento 
Financeiro) que vem sendo implementados no passado recente. 

A situacao mais usual e J c > J a . Para agentes com ampla capacidade de financiamento 
de suas operacoes, o custo de oportunidade para captacao sera considerado igual a taxa 
para aplicacoes. Agentes pouco capitalizados, por outro lado, frequentemente tern que 
recorrer a mecanismos de financiamento para manutencao de suas operacoes na bolsa, 
para pagamento da margem e dos ajustes diarios. Para esses agentes a taxa de captacao 
sera considerada como sendo J c . O problema desses agentes descapitalizados, muitas 
vezes produtores rurais, e um problema serio a ponto de existir interesse pelo 
financiamento para operacoes em bolsa para commodities agrfcolas, considerando taxas 
de juros menores que as praticadas pelo mercado (ex. Funcafe e Fundo Agropecuario BB- 
BM&F). 

Definicao de contas para o fluxo de caixa e obtencao do custo total da operacao - CT 

Para modelar o fluxo de caixa na operacao, consideramos 2 contas, representadas por R e 
A. A primeira e a conta R que contabiliza todos os fluxos monetarios reais da opera§ao, 
sendo atualizada no final de cada perfodo. Nesta conta sao creditados e/ou debitados: (a) 
opera§6es associadas a margem de garantia; (b) taxas operacionais; (c) ajustes diarios; 
(d) tributos; e, (e) encargos financeiros. Quando o saldo e positivo, no dia anterior, essa 
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conta e creditada dos juros obtidos pela aplicacao da taxa J a sobre o saldo anterior. Caso 
o saldo seja negativo, no dia anterior, a conta e debitada de juros a taxa J c sobre esse 
saldo. Para agentes capitalizados, com capacidade propria de financiamento de suas 
operacoes, assume-se J c =J a . 

A segunda e a conta A, que contabiliza somente os ajustes diarios. Essa conta 
possibilitara o calculo do IRPF e CPMF nos saldos dos ajustes (quando devidos), para 
debito na conta R nos momentos apropriados. Em operacoes reais podem existir mais 
contas. As 2 contas consideradas, contudo, caracterizam de uma forma simplificada a 
natureza do processo. A existencia de mais contas pode implicar em custos adicionais, 
nao considerados neste desenvolvimento. As datas correspondentes ao final dos meses, 
quando certos tributos sao calculados, podem ser definidas especificamente ou assumidas 
como ocorrendo a cada 22 dias uteis de operacao (pressuposicao utilizada para obtencao 
dos resultados na secao 5). 



Se representarmos por S m o saldo existente na conta R, ao final do perfodo m, na 
ausencia de custos 

P 0 f =S«+P,l, (4.2) 

ou seja, os ajustes somados ao preco final seriam iguais ao pre§o no perfodo 0, o pre§o 
fixado. Como os custos (CT) existem, a situa§ao mais freqiiente e 

P 0 f =SZ+P£+CI (4.3) 

ou 

CT = (P 0 / -P m / )-^. (4.3.1) 

O termo entre parentesis na expressao (4.3.1) representa exatamente o saldo da conta A, 
observado no perfodo m, logo 

CT = .(4.4) 
Se os custos considerados fossem nulos, os saldos finais das contas A e R seriam iguais. 
Em geral, como esses custos existem, a conta R apresentara, em geral, saldo final menor 
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que o existente na conta A. Uma medida de custo mais facil de ser entendida, pode ser 

f 

obtida pela sua definicao como um percentual de Pq , o preco fixado, por 



Na verdade, CT e uma variavel aleatoria dependente das trajetorias de precos 
consideradas, sendo uma fun5ao dessas trajetorias, podendo ser caracterizada pelo 
processo de simula§ao desenvolvido na secao 3. Esse custo pode ser decomposto, do 
ponto de vista contabil, em seus diferentes componentes: juros nominais recebidos(JR), 
juros nominais pagos (JP), taxas da bolsa (TB), IRPF, CPMF, ajustes (A), compra do 
certificado ( CDB 0 ) e venda do certificado (CDB m ). 

Para um melhor entendimento do processo de decomposi§ao dos custos em seus varios 
componentes, utilizado no trabalho, e apresentada, a seguir, a evolucao do saldo da conta 
R, em cada periodo, pelas equacoes 



^=-CDB 0 -TB 0 

S? =S*+A l + JRj - JPj - IRP^ - CPMIj 

sLi = S«_2 + A m .! + JR m .! - JP m4 - IRPF m _, - CPMF^.j 

Si = Ci + A m - TB m + JR m - JP m - IRPF m - CPMF m + CDB m . (4 .6) 



Para simplificar o modelo, assume-se que o primeiro CPMF e pago no periodo 1 da 
opera§ao e o ultimo no proprio periodo m. A soma das expressoes em (4.6) resulta em 



S£ =-CDB 0 +CDB m -TB 0 -TB m + ^ =i (A t + JR t - JP t -IRPF t -CPMF t ). .(4.7) 



Juntando as expressoes (4.4), (4.7) e (4.8), chega-se a defini§ao do custo CT a partir de 
seus componentes por: 




xlOO. 



.(4.5) 



No caso da conta A, o saldo no periodo m seria definido simplesmente por 




.(4.8) 



50 



CT = CDB 0 - CDB m + TB 0 + TB m - ^ JR t + ^ JP t + ^IRPF t + £^CPMF t ^ 

Para uma melhor caracterizacao dos custos associados a aquisicao do CDB para cobertura 
da margem, consideramos, inicialmente, que o ganho direto, em juros recebidos pela 
operacao com o CDB, e definido por 

JCDB R - CDB m - CDB 0 . (4. 10) 

c 

Por outro lado, o custo do carregamento desse CDB, representado aqui por JCDB , 
dependera, da situacao do agente com relacao a disponibilidade de fundos. Se for um 
agente com ampla disponibilidade de recursos, o custo sera definido pela taxa J a , o custo 
de oportunidade do capital do agente para aplicacoes, pelo penodo m. Como assumimos 
que J a e a propria taxa de juros associada a aplicacoes em CDBs, ocorrera nesse caso 
que 

JCDB R =JCDB C (4.11) 

o que implica a inexistencia de custos para manutencao da margem nesse caso. 
Alternativamente, uma outra situacao extrema seria a de um agente com total falta de 
capital, tendo que captar no mercado esses recursos a taxa J c . Nesse caso, J c > J a , de 
onde decorre que 

JCDB R <JCDB C . (4.12) 

C R 

A diferen§a JCDB -JCDB corresponde ao custo lfquido de carregamento desse CDB. 
Se J ^ e representam taxas diarias, esse custo seria definido por 

JCDB C -JCDB R =CDB 0 (l + / c rf ) m -CDB 0 (l + /f) ra . (4.13) 

Para facilitar a decomposi§ao dos custos de forma a isolar o componente relacionado a 
margem de garantia, fazemos a substitui§ao de (4.10) em (4.9), e somamos e subtramios 

JCDB" na expressao resultante, levando a 
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CT = JCDB C -JCDB C -JCDB* +TB 0 + TB m - JR t + 2Zi JPt + ZTi IRPFt + 2]"i CPMFt 

(4.14) 

Se definirmos, para simplificar ultima expressao, 

CMARGEM= JCDB C - JCDB R , (4.15) 
como o custo da margem, associada ao custo de carregamento do CDB, 

TB = TB 0 + TBj (4.16) 

como o total de taxas operacionais, 

tj =-jcdb c -xr =1 jR t +zr=i jPt ( 4 - 17) 

como o total de juros pagos, deduzidos os juros pagos para manutencao da margem, 

TIRPF= IRPF t (4.18) 

como o total pago de imposto de renda (IRPF), e 

TCPMF = CPMF t (4.19) 

como o total de CPMF pago na operacao. A partir dessas ultimas expressoes, o custo total 
(CT), apresentado em (4.14) seria redefinido por 

CT = CMARGEM+ TB + TJ + TIRPF+ TCPMF. (4.20) 

Nessa expressao, o custo total CT foi decomposto em cinco componentes, que tambem 

podem ser expressos na forma percentual, com relacao a > tomando por base o 
procedimento utilizado na expressao (4.5). 



Definicao dos algoritmos para obten^ao dos componentes de CT 

A especificacao formal dos procedimentos utilizados para defini§ao de cada componente 
em (4.20), sera apresentada a seguir, na forma de algoritmos. As informa§6es requeridas 
para a obten§ao desses componentes, para cada trajetoria de pre§os, incluem: (a) 
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definicao da trajetoria de precos Pq ,P\ ,"',P m (a partir de procedimento descritos na 
Secao 3), considerando que os precos sao os valores de 1 contrato; (b) taxa operacional 
basica, taxa de registro e emolumentos, representados nos algoritmos pelas variaveis 
taxatob, valreg e taxaemol; (c) taxas de juros diarias J d a e J d c ; (d) valores da margem 
de garantia, representado por valormargem (e) alfquotas do IRPF e da CPMF, 
representadas nos algoritmos por taxairpf e taxacpmf. Deve especificar, tambem, o tipo 
de agente com relacao a disponibilidade de fundos, na variavel tipoagente. Para agentes 
que dispoe de recursos amplos para financiamento proprio das operacoes, considera-se 
tipoagente=l, para agentes que precisam captar recursos para operacao, considera-se que 
tipoagente=2. 

Inicialmente, o algoritmo 4A (Figura 4.1), define os saldos diarios na conta A. A partir 
dessas definicao, o algoritmo 4B (Figura 4.2), especifica o calculo dos valores devidos 
de IRPF a cada fim de mes, associados a operacao no mercado futuro. No mundo real o 
IRPF deve pago 15 dias uteis apos o final do mes. No modelo considerado, o pagamento 
e considerado no proprio final do mes com um desconto considerado. 

A seguir, o algoritmo 4C (Figura 4.3), apresenta a definicao dos saldos diarios da conta 
R, assim como dos componentes dos custos da operacao, especificados em (4.20). 
Finalmente, o algoritmo 4D (Figura 4.4) apresenta uma descricao simplificada do 
processo completo de simulacao para um contrato de venda. Para um contrato de compra 
basta alterar o calculo do ajuste (como indicado na figura). O algoritmo considera os 
procedimentos descritos na Secao 3 para obtencao de trajetorias de precos a partir de 
modelos alternatives. Como resultado, apresenta uma amostra tamanho n do vetor que 
caracteriza os componentes do custo da opera§ao, condicionados pelas trajetorias de 
precos simuladas. Os algoritmos apresentados indicam uma possfvel seqiiencia de 
procedimentos para implementa§ao das ideias apresentadas nesta secao e na se§ao 
anterior. Foram desenvolvidos com o interesse de clarificar os procedimentos, e nao com 
o intuito de buscar a maxima eficiencia computacional, o que poderia dificultar seu 
entendimento. 
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Algoritmo 4A - Def inigao dos Saldos da Conta A 

inicio 

especifique: m (numero de dias da operagao) 

P<f , P{ , • ■ (trajetoria de pregos) 

para i<—\,---,m, faga{ 

5, A <- S fi x + (P/ - Pf_ x ) (inverta pregos se for contrato de 

compra) 

} 

intorme jq >"i 



Figura 4.1 Algoritmo utilizado para definicao dos saldos da conta A 



Algoritmo 4B - Calculo do Imposto de Renda 

inicio 

especifique: k (num. de fins de mes na operagao) 

d l ,d 2 ,---,d k (num. dos dias de fim de mes) 

taxairpf (alfquota do imposto) 

Sq , Sj , • • •, (saldo da conta A nos dias de operagao) 

J (taxa de juro para calculo do desconto 

por pagamento antecipado) 

d 0 <— 0 
A^O 

para /<— \,---,k , faga{ 

SM; <— s£ - Sj^ (definigao dos saldos mensais)) 

} 

para i <— l,---,k , faga { 
A^A + SM; 

irpf ,.<-0 
se A > 0 { 

irpf ( A X taxairpf) /(l + J ) 1 5 

A^O 

} 

} 

informe imposto i ,i = l,---,k (total e valores mensais de irpf) 

fim 

obs.: procedimento desenvolvido segundo IN da SRF 25/01 



Figura 4.2 Algoritmo utilizado para compute do imposto de renda (IRPF) a pag; 



Algoritmo 4C- Definigao dos saldos da conta R e de CT 

inicio 

especifique: m, Pq , P/ , ■ • ■, Pj, (dias de oper. e trajetoria de pregos) 

d l ,d 2 ,---,d k (num. dos dias de fim de mes) 

irpf j, j = l,---,k (valores mensais de irpf) 

tipoagente (tipo do agente quanto a liquidez, 1 ou 2) 

J d , J d (juros para captagao e aplicagao) 

valormargem (valor da margem para o contrato) 
taxatob, valreg, taxaemol (fob, registro, emolumentos) 

taxacpmf, S m (aliquota da cpmf, saldo conta A em m) 

TB 0 <— pj x taxatobx (1 + taxaemol) + taxareg (taxas bolsa) 
CDB 0 <— valormargem 

Sq ^-TB 0 -CDB 0 
CPMF !<- TB 0 X taxacpmf 
JP,JR,TIRPF,j^O 

se tipoagente=2 faga { J <— J d c } caso contrario faca { J <— J d } 
para i <- \,---,m , faga{ 

se S£j > 0 faga { juros <- s£j x J d a ; JR <- JR + juros } 
caso contrario faga { juros <— S*j x J ■ JP <— JP - y'wros } 

Sf <- Sf_ x + (fj/j - P/ ) +juros (inverta pregos em contr. compra) 
se z'e {d x ,d 2 ,---,d k } faga{ j<—j + \ 

Sf -irpfj 
TIRPF <— TIRPF+ irpf ■ j 

} 

} 

TB m x taxatobx (1 + taxaemol) + taxareg (taxas bolsa) 

CDB m <- valormargemx(l + 7^) ffl 
CPMF m <- TB m X taxacpmf 

seS*<0 faga { CPMF m <— CPMF m + taxacpmf x (S^ ) } 

S^S*-TB m +CDB m -CPMF m 

TB<— TB 0 +TB m 

CMARGEM<— valormargemx((l + 7) m -{l + J d a ) m ) (custo da margem) 
TCPMF <— CPMF 0 +CPMF m 

TJ <— JP - JR - (valormargemx (1 + J) m - valormargem) 

CT <- CMARGEM+ TB + TJ + TIRPF+ TCPMF 
informe CT, CMARGEM, TB, TJ, TIRPF, TCPMF 

fim 

Figura 4.3 Algoritmo utilizado para compute do saldo na conta R 
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Algoritmo 4D - Processo Completo (Simplificado) 

inicio 

especifique: n (num. de simulagoes) 

m, (dias de oper. e trajetoria de pregos) 

d l ,d 2 ,---,d k (num. dos dias de fim de mes) 

tipoagente (tipo do agente quanto a liquidez, 1 ou 2) 

J d , J d (juros para captagao e aplicagao) 

tipoagente (tipo do agente quanto a liquidez, 1 ou 2) 
valormargem (valor da margem para o contrato) 
taxatob, valreg, taxaemol (tob, registro, emolumentos) 
taxacpmf (alfquota da cpmf, saldo conta A em m) 
M (modelo do processo de simulagao de pregos) 

P , 9 (prego em 0 e parametros para o modelo M) 

para s <- l,---,n , faga{ 

P({ , P/ , ■ ■ ■, <— simule a trajetoria de pregos s com M, P , 9 

Sq , 5j A , • • • , <— obtenha saldos da conta A com Algoritmo 4 A 
irpf\ • • -irpf k <— obtenha o irpf no final dos meses pelo Algoritmo 4B 

CMARGEM^B^^TIRPI^TCPMF 6 <- obtenha os componentes 

do custo para a trajetoria s pelo algoritmo 4C 
CMARGEM^ , TB \ , TJ \ , TIRPE^ , TCPME^ <- 

CMARGEM 8 ,TB S ,TJ S , TIRPF 5 , TCPMF 5 x -i- x 1 00 
(custos definidos na forma percentual) 

CT4 <- CMARGEM^ +TB S % +TJ S % + TIRPF^ + TCPME^ 
pf-pf 

VP% < ^— j — xlOO (variagao percentual do prego entre 0 em) 

} 

obtenha, para as variaveis de interesse, media, intervalo de confianca e 
matriz de correlacao (amostral) a partir de 

CT^CMARGEM^^TB^^TJ^^TIRPE^^TCPME^^V/^, s = \,---,n 

fim 



Figura 4.4 Algoritmo utilizado para computo dos custos operacionais e seus componentes 
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5. Estimativa de custos operacionais em mercados futuros 

Esta secao ilustra a aplicacao dos procedimentos desenvolvidos nas secoes anteriores. 
Um estudo de caso e argumentos teoricos sao utilizados para o entendimento da 
magnitude e natureza de certos custos operacionais em mercados futuros. 

O estudo de caso considera contratos do mercado futuro de milho, transacionados na 
BM&F. Inicialmente, os custos sao calculados a partir da propria evolucao dos precos 
observada no mercado e dos conceitos e algoritmos definidos na secao anterior. Num 
segundo momento, sao estimados modelos que tentam caracterizar, do ponto de vista 
estatfstico, as trajetorias que esses precos podem assumir. A partir dessa caracterizacao, 
realizada com grau crescente de complexidade, sao obtidas estimativas dos componentes 
dos custos, com o uso da tecnica de simulacao, detalhada na Secao 3. 

Para efeito do exercfcio realizado, considera-se que a taxa de juros por dia de operacao 
para aplicacoes (a mesma considerada para o CDB), representada por J d , e 0,00055 
(aproximadamente 1,2% ao mes, considerando 22 dias de operacao por mes). Assume-se 
que a taxa de juros diaria de captacao, representada por J d , e 0,0017 (aproximadamente 
3,75% ao mes). Consideram-se 2 tipos extremos de investidor, o primeiro, tipo 1, 
denominado agente capitalizado, tern plena capacidade de financiar a operacao com 
recursos proprios, considerando um custo de oportunidade igual a J d a . Por outro lado, o 
agente tipo 2, denominado agente descapitalizado, necessita de capital externo para 
custear as operacoes, ao custo de J d c , e aplica os saldos positivos a taxa J d . Considera- 
se tambem uma restricao de oscila§ao maxima diaria dos pre§os (usada pela BM&F), 
definida de forma que o valor absoluto da taxa de varia§ao nao pode exceder 5%. 

Contratos de Milho - BM&F - Custos de trajetorias observadas 

Sao considerados na analise os contratos futuros de milho com vencimento em MAI02, 
JAN03, MAR03 e MAI03, por serem os mais proximos do perfodo de colheita do 
produto. Desses contratos, foram aproveitados os ultimos 121 pre§os de ajuste 
observados, desprezando-se as informa§6es dos primeiros meses do contrato, que 
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usualmente estao associadas a um baixo volume de negociacao. Contratos de milho 
anteriores a MAI02 tambem foram desprezados por 2 razoes principals: (a) 
representavam um outro modelo de contrato; ou (b) o volume negociado atingia mveis 
pouco expressivos. Por outro lado, contratos com vencimento posterior a MAI03 foram 
desconsiderados porque ainda estavam distantes do vencimento na data da pesquisa. A 
Figura 5.1 mostra a evolucao da trajetoria dos precos observada durante o periodo de 
analise considerado para cada um desses contratos. 
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Figura 5.1 Evolucao dos pre§os em contratos futuros de milho selecionados 
nos ultimos 121 dias do contrato (veja texto). 
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Os custos foram definidos de forma percentual (com relacao ao valor do contrato no 
periodo 0) e decompostos por categoria, para cada um desses vencimentos, 
considerando-se um periodo de 120 dias. Esses custos, para investidores do tipo 1 e 2 
(capitalizados e descapitalizados) sao apresentados nas Tabelas 5.1 e 5.2, 
respectivamente. 



Componente 


Custos como percentagem do valor inicial do contrato (120 dias) 


MAI02 


JAN03 


MAR03 


MAI03 


Custo de Margem 


0,00% 


0,00% 


0,00% 


0,00% 


Taxas da Bolsa 


0,71% 


0,83% 


0,69% 


0,63% 


Juros lfquidos 


0,09% 


3,12% 


-0,03% 


-0,51% 


IRPF 


1,75% 


0,00% 


1,82% 


3,11% 


CPMF 


0,08% 


0,23% 


0,06% 


0,002% 


Custo Total 


2,63% 


4,18% 


2,55% 


3,23% 



Tabela 5.1 Custos operacionais estimados em trajetorias de precos observadas no 
contrato de milho (BM&F) para agentes tipo 1 (capitalizado). 



Componente 


Custos como percentagem do valor inicial do contrato (120 dias) 


MAI02 


JAN03 


MAR03 


MAI03 


Custo de Margem 


1,72% 


1,27% 


0,99% 


0,92% 


Taxas da Bolsa 


0,71% 


0,83% 


0,69% 


0,63% 


Juros lfquidos 


0,25% 


9,86% 


-0,05% 


-1,17% 


IRPF 


1,75% 


0,00% 


1,82% 


3,11% 


CPMF 


0,08% 


0,23% 


0,06% 


0,002% 


Custo Total 


4,51% 


12,20% 


3,52% 


3,49% 



Tabela 5.2 Custos operacionais estimados em trajetorias de precos observadas no 
contrato de milho (BM&F) para agentes tipo 2 (descapitalizado). 
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As Tabelas 5.3 e 5.4 apresentam os mesmos custos para uma operacao de 60 dias, 
contados do final do contrato. Essa situa§ao e apresentada para evidenciar a 
dependencia de cada componente do custo no prazo da operacao, com base nos dados 
empfricos de trajetorias. Os itens dependentes de taxas de juros sao claramente mais 
afetados que outros. De um modo geral os custos percentuais sao menores que os 
observados na operacao de 120 dias. Conclusoes mais definitivas, porem, dependerao 
de outras analises que serao apresentadas posteriormente. 



Componente 


Custos como percentagem do valor inicial do contrato (60 dias) 


MAI02 


JAN03 


MAR03 


MAI03 


Custo de Margem 


0,00% 


0,00% 


0,00% 


0,00% 


Taxas da Bolsa 


0,72% 


0,64% 


0,74% 


0,67% 


Juros lfquidos 


0,22% 


0,14% 


0,23% 


0,29% 


IRPF 


0,00% 


1,56% 


0,21% 


0,00% 


CPMF 


0,09% 


0,01% 


0,11% 


0,04% 


Custo Total 


1,03% 


2,35% 


1,28% 


1,00% 



Tabela 5.3 Custos operacionais estimados em trajetorias de precos observadas no 
contrato de milho (BM&F) para agentes tipo 1 (capitalizado). 



Componente 


Custos como percentagem do valor inicial do contrato (60 dias) 


MAI02 


JAN03 


MAR03 


MAI03 


Custo de Margem 


0,82% 


0,38% 


0,51% 


0,50% 


Taxas da Bolsa 


0,72% 


0,64% 


0,74% 


0,67% 


Juros lfquidos 


0,68% 


0,47% 


0,67% 


0,90% 


IRPF 


0,00% 


1,56% 


0,21% 


0,00% 


CPMF 


0,09% 


0,01% 


0,11% 


0,04% 


Custo Total 


2,31% 


3,06% 


2,23% 


2,10% 



Tabela 5.4 Custos operacionais estimados em trajetorias de precos observadas no 
contrato de milho (BM&F) para agentes tipo 2 (descapitalizado). 
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Contratos de Milho (BM&F) - Caracteriza^ao estatistica das trajetorias 

Nesta parte do trabalho as series de precos sao caracterizadas por modelos estatfsticos. 
Esses modelos podem oferecer uma representacao relativamente simples, mas 
razoavelmente adequada, do comportamento estocastico dos precos. A analise 
desenvolvida possibilita conclusoes mais gerais sobre a natureza dos custos em operacoes 
em mercados futuros, que aquelas estabelecidas por uma unica trajetoria. Essa trajetoria 
observada seria uma das infinitas trajetorias consideradas a partir modelo estimado. 



O objetivo central dessa analise e a caracterizacao de um modelo ARIMA(c,i,d)- 
GARCH(p,q) que possa representar satisfatoriamente a evolucao das taxa diaria de 
crescimento dos precos r t , obtida por (2.2) e (2.3). Na Figura 5.2 sao apresentados 
graficos que mostram as evolucoes dessas taxas, para os contratos de milho descritos 
anteriormente. 
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Figura 5.2 Evolucao das taxas de varia§ao dos pre§os nos contratos futuros 

de milho selecionados nos ultimos 120 dias do contrato (veja texto). 
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Os graficos apresentados foram construfdos com a mesma escala, de forma a facilitar a 
comparacao da evolucao das taxas nos quatro contratos. Para facilitar o diagnostico de 
problemas de autocorrelacao (AR) e dependencias em medias moveis (MA), sao 
apresentadas, nas Figuras 5.3 e 5.4, respectivamente, as funcoes de autocorrelacao e 
autocorrelacao parcial das taxas em cada uma das series. Nas figuras, barras que excedem 
as linhas tracejadas indicam rejeicao da hipotese da nulidade para a autocorrelacao ao 
nfvel de 5%. Para as series JAN03, MAR03 e MAI03, indfcios de possfveis efeitos 
associados a AR e/ou MA parecem mais evidentes (veja Mills, 1997 e Hamilton, 1994, 
para detalhes sobre o diagnostico). Adicionalmente, a Figura 5.2 indica que as series 
apresentam periodos (clusters) com maior ou menor volatilidade, e sugere que os 
modelos GARCH podem ser apropriados para modelar o comportamento da variancia 
condicional das taxas. Esses indfcios sao reforcados pela funcao que apresenta a 
autocorrelacao das taxas elevadas ao quadrado, na Figura 5.5 (Andersen et al., 2002 e 
Laurent e Peters, 2002). 
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Figura 5.3 Funcoes de autocorrela§ao das taxas diarias, obtidas das series de pre§os 
de milho selecionados nos ultimos 121 dias do contrato (veja texto). 
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Figura 5.4 Funcoes de autocorrelacao partial das taxas diarias, obtidas 
das series de precos de milho selecionados nos ultimos 121 dias 
do contrato (veja texto). 
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Figura 5.5 Fun§oes de autocorrelacao das taxas diarias ao quadrado, obtidas 
das series de precos de milho selecionados nos ultimos 121 dias 
do contrato (veja texto). 
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Foram testados, para cada serie, modelos ARIMA(c,i,d) que consideraram todas as 
combinacoes possfveis considerando c variando de 0 a 3, i variando de 0 a 1 e d variando 
de 0 a 3. Os modelos testados foram classificados atraves do criterio BIC (Bayes 
Information Criterium), desenvolvido por Schwarz (1978). Esse criterio tende a ser 
considerado mais adequado do ponto de vista teorico - privilegiando modelos com menos 
parametros, por exemplo - que outros criterios populares tambem fundamentados no 
valor da funcao de verossimilhanca como o AIC (Akaike Information Criterium), 
argumenta Mills (1997). A Tabela 5.5 apresenta os 3 melhores modelos sob esse criterio 
(quanto menor o valor do BIC, melhor e o modelo). 



Serie 


Modelos 


Criterio 


de taxas 


Melhores 


BIC 




ARIMA(0,0,0) 


-747,09 


MAI02 


ARIMA( 1,0,0) 


-742,66 




ARIMA(0,0,1) 


-742,60 




ARIMA( 1,0,0) 


-648,27 


JAN03 


ARIMA(0,0,1) 


-646,44 




ARIMA( 1,0,1) 


-643,22 




ARIMA(2,0,0) 


-660,58 


MAR03 


ARIMA(0,0,2) 


-660,08 




ARIMA(0,0,1) 


-659,98 




ARIMA(2,0,0) 


-710,30 


MAI03 


ARIMA(0,0,1) 


-710,18 




ARIMA(3,0,0) 


-707,10 



Tabela 5.5 Melhores modelos ARIMA para representacao de cada serie, 
segundo o criterio BIC (Schwarz, 1978) - veja texto. 
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Os resfduos do melhor modelo ARIMA obtido em cada serie foram utilizados para a 
selecao de um modelo GARCH(p,g). Os modelos testados consideraram todas as 
possibilidades com p variando de 0 a 2 e q variando de 0 a 2. Os modelos testados foram 
classificados atraves do criterio BIC e estao apresentados em ordem de superioridade por 
esse criterio na Tabela 5.6. 



Serie 


Modelos 


Criterio 


de taxas 


Melhores 


BIC 




GARCHfO 1) 


-751,63 


MAI02 


GARCH(1,1) 


-747 ,47 




GARCH(0,0) 


-747,10 




GARCH(1,1) 


-654,61 


JAN03 


GARCH(0,1) 


-653,43 




GARCH(0,0) 


-653,22 




GARCH(0,1) 


-750,94 


MAR03 


GARCH(1,1) 


-744,92 




GARCH(0,2) 


-738,11 




GARCH(0,1) 


-678,00 


MAI03 


GARCH(1,1) 


-677,00 




GARCH(0,2) 


-671,84 



Tabela 5.6 Melhores modelos GARCH para representa§ao de cada serie, 
segundo o criterio BIC (Schwarz, 1978) - veja texto. 
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Os modelos ARIMA e GARCH mais promissores, escolhidos a partir da analise 
preliminar realizada nos ultimos paragrafos, foram testados em conjunto, tambem pelo 
criterio BIC. Dessa analise foram definidos, na Tabela 5.7, modelos adequados para 
representar as series de precos consideradas. Nessa Tabela sao especificados os 
parametros estimados pelo metodo da maxima verossimilhanca, com o apoio do Software 
R, versao 1.7.1 e do software Ox-G@rch 3.0. Nos modelos especificados para MAI02 e 
MAI03, a distribuicao t-Student, respectivamente com 10 e 6 graus de liberdade, 
corrigida para apresentar uma variancia unitaria, foi utilizada em lugar da normal 
padronizada, para caracterizar a distribuicao de z em (2.11). Essa pressuposicao foi 
suportada pelo criterio BIC, em compara§ao a outras pressuposicoes (ex. normalidade). 
Nos outros casos, a distribui§ao normal padronizada foi utilizada. A especifica§ao dos 
parametros nos modelos ARIMA e GARCH segue a nota§ao utilizada em (2.7), (2.8) e 
(2.11). As Figuras Al a A4, no Apendice, ilustram algumas trajetorias de taxas 
amostradas desses modelos, obtidas a partir do algoritmo 3D, apresentado na secao 3. A 
simula§ao das trajetorias utilizou os parametros estimados como sendo os verdadeiros 
valores dos parametros, uma pratica comum em estudos desse tipo - veja por exemplo 
Laurent e Peters (2002). Um procedimento mais elaborado e tecnicamente mais correto, 
porem menos utilizado, poderia considerar toda a informacao contida na distribui§ao a 
posteriori dos parametros (Lence, 1995). Para cada trajetoria, seriam simulados os 
valores dos parametros dessa distribui§ao. Esse procedimento, contudo, nao e utilizado 
neste estudo. 
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Serie 
de taxas 


Estimativas 
dos Parametros 


Modelo 


IVIAIUZ 


a 0 =0,00008124 
b x = 0,3839 


ARIMA(0,0,0) 


JAN03 


8j =0,4614 
a 0 =0,00003164 
a x =0,5510 

b x =0,3415 


ARIMA( 1,0,0) 
GARCH(1,1) 


MAR03 


§! =0,4514 

§ 2 = -0,1401 
a 0 = 0,00007539 
a x = 0,4234 

b x = 0,2785 


ARIMA(2,0,0) 
GARCH(1,1) 


MAI03 


§! =0,3688 
9 2 = -0,05567 
a 0 = 0,00004488 
a, = 0,2508 

fci =0,5513 


ARIMA(2,0,0) 
GARCH(1,1) 



Tabela 5.7 Modelos selecionados para representacao de cada serie. 



Contratos de Milho (BM&F) - Estimativa dos Custos pela Simula^ao 

As Tabelas 5.8 e 5.9 apresentam as medias dos componentes dos custos, respectivamente, 
para o agente tipo 1 (capitalizado) e tipo 2 (descapitalizado). Essas medias foram obtidas 
a partir de 10.000 trajetorias de precos simuladas atraves dos procedimentos descritos nas 
secoes anteriores. Os resultados tendem a apresentar um quadro com menos extremos, em 
funcao de que os custos representam a media das trajetorias possfveis obtidas dos 
modelos estimados. Nas trajetorias efetivamente observadas, analisadas anteriormente, os 
precos, em geral, apresentaram tendencia de crescimento, o que aumenta a necessidade 
de desembolsos para os ajustes e reduz a magnitude do imposto de renda devido (em 
funcao dos saldos negativos na conta A). Na simula§ao esse problema e minimizado pela 
consideracao de uma amostra de trajetorias mais ampla, a partir de um modelo estatfstico 
selecionado para representar a natureza fundamental das trajetorias. Para simplificar o 
modelo, assume-se nesses experimentos que a taxa de desconto para pagamento 
antecipado em 15 dias do IRPF e zero. 4 

A incerteza existente com relacao a trajetoria de precos acaba influenciando a propria 
distribuiijao de custos. Na Figura 5.6 sao apresentadas estimativas da distribui§ao de 
probabilidade do custo, para os modelos correspondentes a MAI02 e JAN03, para o 
agente tipo 1 (capitalizado). A Figura 5.7 apresenta estimativas da distribui§ao de 
probabilidade dos custos para os mesmos vencimentos para agentes tipo 2 
(descapitalizados). O objetivo dessas figuras e o de ilustrar a natureza da incerteza 
existente. 

As Tabelas 5.10 e 5.11 contem as matrizes de correlacao simples obtidas a partir das 
10.000 trajetorias simuladas. Essas matrizes descrevem a associa§ao estatfstica entre os 
diferentes componentes do custo, para as trajetorias simuladas pelos modelos 
correspondentes a MAI02 e JAN03, para o agente tipo 2 (capitalizado). 



4 O impacto dessa pressuposicao e mmimo, conforme indicaram experimentos preliminares: aumento, em 
menos de 1% em termos relativos, da magnitude do custo total medio em operacoes de 120 dias, quando 
comparado a situacoes em que o desconto e considerado. 
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Para verificar a sensibilidade dos resultados ao refinamento do modelo, e realizado um 
novo experimento, que considera um modelo extremamente simples para caracterizacao 
de cada serie, correspondente ao caso ARIMA(0,0,0)-GARCH(0,0), que se reduz ao 
modelo descrito em (2.2) a (2.4). A volatilidade diaria foi estimada, nesse caso, pelo 
simples desvio padrao das taxas observadas no periodo. Essas estimativas obtidas, para o 
agente tipo 2 sao apresentadas na Tabela 5.12. Os custos obtidos pelo modelo 
simplificado (Tabela 5.13) sao menores que os obtidos a partir da utilizacao do modelo 
mais elaborado, o que sugere que modelos mais aderentes a realidade do comportamento 
observado dos precos, podem contribuir para estimativas mais apropriadas (e 
eventualmente maiores) dos custos operacionais. Deve-se ressaltar, entretanto, que apesar 
das limitacoes, esse modelo simplificado pode proporcionar aproximacoes 
suficientemente boas para certas aplicacoes. 
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Componente 


Custos como percentagem do valor inicial do contrato (120 dias) 


MAI02 


JAN03 


MAR03 


MAI03 


Custo de Margem 


0,00% 


0,00% 


0,00% 


0,00% 


Taxas da Bolsa 


0,64% 


0,63% 


0,65% 


0,64% 


Juros lfquidos 


0,07% 


0,13% 


0,20% 


0,03% 


IRPF 


1,41% 


3,67% 


2,63% 


2,35% 


CPMF 


0,02% 


0,04% 


0,04% 


0,03% 


Custo Total 


2,14% 


4,47% 


3,52% 


3,05% 



Tabela 5.8 Custos operacionais medios estimados em 10.000 trajetorias de 
precos simuladas a partir de modelos estimados para agentes 
tipo 1 (capitalizados). 



Componente 


Custos como percentagem do valor inicial do contrato (120 dias) 


MAI02 


JAN03 


MAR03 


MAI03 


Custo de Margem 


1,72% 


1,27% 


0,99% 


0,92% 


Taxas da Bolsa 


0,64% 


0,63% 


0,64% 


0,64% 


Juros lfquidos 


0,19% 


0,68% 


0,79% 


0,47% 


IRPF 


1,47% 


3,86% 


2,74% 


2,36% 


CPMF 


0,02% 


0,04% 


0,04% 


0,03% 


Custo Total 


4,04% 


6,48% 


5,21% 


4,41% 



Tabela 5.9 Custos operacionais medios em 10.000 trajetorias de precos simuladas 
a partir de modelos estimados para agentes tipo 2 (descapitalizados). 
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CT 


CMARGEM 


TB 


TJ 


TIRPF 


TCPMF 


Varia^ao 


CT 


1,00 














CMARGEM 


0,00 


1,00 












TB 


0,13 


0,00 


1,00 










TJ 


0,39 


0,00 


0,85 


1,00 








TIRPF 


0,19 


0,00 


-0,84 


-0,83 


1,00 






TCPMF 


0,47 


0,00 


0,87 


0,76 


-0,54 


1,00 




Varia^ao 


0,13 


0,00 


1,00 


0,84 


-0,84 


0,87 


1,00 



Tabela 5.10 Matriz de correlacao amostral entre os componentes do custo do 
modelo associado contrato MAI03 estimados em 10.000 
trajetorias de preco simuladas para agentes tipo 2 (descapitalizados). 





CT 


CMARGEM 


TB 


TJ 


TIRPF 


TCPMF 


Varia^ao 


CT 


1,00 














CMARGEM 


0,00 


1,00 












TB 


0,21 


0,00 


1,00 










TJ 


0,49 


0,00 


0,85 


1,00 








TIRPF 


0,21 


0,00 


-0,81 


-0,75 


1,00 






TCPMF 


0,54 


0,00 


0,88 


0,78 


-0,48 


1,00 




Varia^ao 


0,21 


0,00 


1,00 


0,85 


-0,81 


0,88 


1,00 



Tabela 5.11 Matriz de correlacao amostral entre os componentes do custo do 
modelo associado contrato JAN03 estimados em 10.000 
trajetorias de preco simuladas para agentes tipo 2 (descapitalizados). 
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Figura 5.6 Distribui§6es de probabilidade de custos CT (%) estimadas, a partir 
das trajetorias simuladas para MAI02 e JAN03, para o agente tipo 1 
obtidas pelo modelo simplificado. 
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Figura 5.7 Distribuicoes de probabilidade de custos CT (%) estimadas, a partir 
das trajetorias simuladas para MAI02 e JAN03, para o agente tipo 2 
obtidas pelo modelo simplificado. 
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Serie 


Volatilidade 
diaria 


Volatilidade 

Diaria 
Anualizada 


MAI02 


0,0109 


0,1737 


JAN03 


0,0172 


0,2741 


MAR03 


0,0176 


0,2804 


MAI02 


0,0135 


0,2152 



Tabela 5.12 Volatilidade dos contratos estimada atraves da trajetoria 
real dos precos (veja texto) 



Componente 


Custos como percentagem do valor inicial do contrato (120 dias) 


MAI02 


JAN03 


MAR03 


MAI03 


Custo de Margem 


1,72% 


1,27% 


0,99% 


0,92% 


Taxas da Bolsa 


0,64% 


0,65% 


0,63% 


0,64% 


Juros lfquidos 


0,20% 


0,53% 


0,52% 


0,43% 


IRPF 


1,40% 


2,15% 


2,31% 


1,76% 


CPMF 


0,02% 


0,03% 


0,03% 


0,03% 


Custo Total 


3,99% 


4,65% 


4,50% 


3,77% 



Tabela 5.13 Custos operacionais estimados em 10.000 trajetorias de precos simuladas 
a partir do modelo simplificado, para agentes tipo 2 (descapitalizados). 
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Efeito da volatilidade e numero de dias no custo da operacao 

Os resultados apresentados mostram que o numero de dias da operacao e a volatilidade 
podem afetar a magnitude e a incerteza associadas ao custo de uma operacao no mercado 
future Os efeitos dessa dependencia sao explorados a seguir, a partir de um modelo 
teorico que considera que as taxas de variacao dos precos seguem o modelo basico 
especificado em (2.2) a (2.3), correspondente ao caso simples ARIMA(0,0,0)- 
GARCH(0,0). Nesse modelo, as trajetorias sao totalmente determinadas pela volatilidade 
e pelo valor inicial do contrato. Esse modelo, explorado anteriormente, pode ser util para 
o entendimento de certos efeitos, a custa de algumas simplificacoes. 

Para os experimentos realizados, a margem de garantia e definida como 6,2% do valor do 
contrato, as taxas de juros (a.m.) definidas em J a = 1,2% e J c = 3,75% , a TOB em 0,3%, os 
emolumentos em 6,32% da TOB e a taxa de registro como sendo R$ 0,20 por contrato. 
Foram considerados nos experimentos volatilidades (anualizadas) variando de 0,15 a 0,35 
(15 % a 35%) e operacoes variando de 10 a 120 dias. Os custos totais medios, obtidos a 
partir de 30 mil simulacoes, expressos como um percentual do valor inicial do contrato, 
para os agentes tipo 1 (capitalizado) e tipo 2 (descapitalizado), sao apresentados na 
Figura 5.8. A Figura 5.9 apresenta a estimativa da distribuicao de probabilidade do custo 
total para uma volatilidade anualizada de 25%, em operacoes com duracoes distintas. Na 
Figura 5.10, sao apresentados os custos medios de IRPF, expressos como um percentual 
do valor inicial do contrato. 

Esses resultados indicam que o custo tende a ser mais elevado em situacoes de maior 
volatilidade e em operacoes mais longas. Um outro aspecto importante e a grande 
incerteza que existe no inicio de uma operacao quanto ao seu custo, mesmo que sejam 
perfeitamente conhecidos os parametros estatfsticos que caracterizam as trajetorias de 
precos. Isso e evidenciado na Figura 5.9. Essa incerteza aumenta com volatilidades 
maiores e operacoes mais longas, se traduzindo num risco de custo , que ao lado do risco 
de base, pode dificultar operacoes pelo impedimento de um hedge que possa 
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efetivamente fixar o preco. Esse problema seria menos serio se custo fosse fixo, 
conhecido exatamente no infcio da operacao, algo que em geral nao ocorre. 




dias 




dias 



Figura 5.8 Custos totais medios como % do valor inicial do contrato 

estimados a partir do modelo basico para agentes tipo 1 e 2 
(cescapitalizado e descapitalizado). 
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Figura 5.9 Estimativas da distribuicao do custo total como % do valor inicial do 
contrato a partir do modelo basico para agente tipo 2 \(capitalizado) 
para uma volatilidade anualizada de 0,25, em funcao do numero de dias 
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Figura 5.10 Custos de IRPF medios como % do valor inicial do contrato, 

estimados para diferentes volatilidades e duracoes de operacao. 
(para agentes tipo 1 e tipo 2) 



Estimativa dos custos de novas operates 



Os procedimentos metodologicos desenvolvidos nas secoes anteriores podem ser 
estendidos para possibilitar estimativa do custo de uma nova operacao. Nesse caso, os 
dados existentes ate o momento da operacao podem ser analisados e a partir deles 
estimados os modelos que caracterizem o processo estocastico associado as taxas de 
varia§ao, com o auxflio, por exemplo, do modelo geral ARIMA-GARCH. As trajetorias 
seriam simuladas pelo prazo previsto para a operacao, considerando as informa§oes 
existentes no imcio da operacao. Estimativas das distribui§oes dos componentes dos 
custos e seus momentos amostrais podem ser obtidos pelos procedimentos anteriormente 
descritos. A discussao de procedimentos especfficos sobre esse processo nao e objeto 
desta pesquisa. 
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6. Considera^oes finais 



O trabalho apresentou procedimentos para caracterizacao e simula<jao de trajetorias de 
precos, fluxos de caixa e custos operacionais, aplicados a operacoes em mercados 
futuros. Os modelos utilizados capturam aspectos relevantes de operacoes reais, com 
algumas simplificacoes que tern o objetivo de facilitar a exposicao e implementacao. 
Esses modelos foram comparados, verificando-se diferencas obtidas com graus distintos 
de complexidade e aderencia a realidade. Os metodos utilizados foram formalizados com 
algoritmos apresentados no trabalho, e implementados em software, visando sua 
aplica§ao efetiva. 

Os custos operacionais estimados neste trabalho referem-se ao caso em que o agente e 
uma pessoa fisica , sendo obtidos a partir de procedimentos que consideram a evolucao do 
fluxo de caixa da operacao, condicionados por trajetorias de precos: (a) observadas em 
contratos futuros, (b) simuladas a partir de modelos estatfsticos (ARIMA-GARCH) que 
capturam aspectos relevantes de trajetorias observadas no mercado e (c) simuladas por 
modelos teoricos que dependem so da volatilidade e do numero de dias da operacao. 
Deve-se destacar que os custos operacionais para outros agentes (ex. pessoa juridica com 
lucro presumido, pessoa juridica com lucro calculado, institui§oes isentas de certos 
impostos, etc.) podem ser diferentes dos custos operacionais apresentados neste trabalho, 
sendo assunto de pesquisas futuras. 

A metodologia desenvolvida foi utilizada para estimativa de custos associados a margem 
de garantia, encargos financeiros, imposto de renda, CPMF e taxas operacionais da bolsa, 
para um agente economico do tipo pessoa fisica, operando no mercado futuro. Uma 
dificuldade existente no calculo desses custos e o fato de que, em larga medida, alguns de 
seus componentes - como o imposto de renda e os encargos financeiros - dependem 
fortemente da trajetaria de precos observada em uma operacao. Outro problema e a 
dependencia desses custos a situacao do agente economico quanta ao seu grau de 
necessidade de capital externo. Para agentes descapitalizados, dependentes da capta§ao 
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de recursos para fazer frente a margens, ajustes, impostos, e gerenciamento dos contratos 
de um modo geral, os custos tendem a ser mais elevados que para agentes que dispoem 
de ampla capacidade propria de recursos financeiros. 

Os pressupostos utilizados na aplica§ao dos impostos abstraem certos aspectos da 
realidade. O relaxamento desses pressupostos, contudo, examinado na pesquisa, mostra 
impacto mrnimo sobre os resultados apresentados. De um modo geral, o impacto maior e 
no sentido da subestimativa da tributacao, especialmente do custo associado a CPMF. 
Operacoes reais podem envolver muitas transferencias entre contas, o que pode resultar 
em tributacao por CPMF maior que a obtida neste estudo. 

Os resultados obtidos sugerem que os custos de operacoes no mercado futuro podem ser 
expressivos. Para uma pessoa fisica, numa operacao de 60 dias, sob pressuposicoes 
minimalistas, os custos operacionais medios atingiram valores entre 1,5 e 2% do valor 
inicial do contrato (Figura 5.8, pag. 75), quando a volatilidade anualizada e 20%. Para 
uma operacao com 120 dias, esse custo operacional medio variou entre 2,5 a 3,5%. Para 
uma volatilidade de 30% os custos medios observados estiveram entre 2,5 a 3% e entre 3 
a 4,5%, para 60 e 120 dias, respectivamente. Os valores em cada caso dependem da 
situa§ao da pessoa fisica com relacao necessidade de capta§ao de recursos - os maiores 
valores sao para agentes descapitalizados. Para obtencao desses resultados foram 
pressupostas taxas de juros nas aplica§oes e nas capta§6es, respectivamente, de 1,2% e 
3,75% ao mes. Em situacoes praticas as taxas de capta§ao podem ser muito superiores a 
esse valor, nas condi§6es atuais do mercado. Do lado das aplica§6es, assume-se que os 
saldos positivos de ajustes podem ser investidos a taxa de aplica§ao considerada (algo 
facilitado pelo FIF). O uso de pressuposicoes sobre a dinamica da evolucao dos precos 
mais proximas a realidade (com modelos mais complexos) levou a estimativas de custos 
operacionais superiores as apresentadas acima, obtidas pelos modelos simplificados. 
Esses modelos simplificados foram, contudo, uteis para caracterizar a natureza do 
impacto da volatilidade e tempo decorrido na operacao nos custos da operacao. Deve-se 
destacar que os custos para outros agentes (pessoa juridica com lucro presumido, pessoa 
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juridica com lucro calculado, instituicoes isentas de certos impostos, etc.) podem ser 
diferentes desses custos apresentados. 

E importante perceber que os custos apresentados acima sao aplicaveis as duas 
contrapartes de uma mesma operacao, tanto na posicao vendida como na posicao 
comprada. Isso significa que o custo da transacao sera o dobro dos valores medios 
obtidos acima. E possfvel ainda que numa operacao os custos operacionais para a posicao 
vendida e a posicao comprada sejam diferentes (assimetricos) em funcao da forma com 
que os encargos financeiros e a tributacao sao aplicados. Numa posicao de venda, uma 
trajetoria de precos ascendente pode motivar o pagamento de um elevado volume de 
encargos financeiros pelos recursos necessarios para que o agente faca frente aos ajustes 
necessarios. Para o contrato de compra, por outro lado, haveria recebimentos de juros 
decorrentes da aplicacao dos ajustes positivos, mas os desembolsos com tributacao 
seriam elevados. Um exemplo extremo, que destaca esses efeitos, analisado neste 
trabalho, e ilustrado pelo contrato do milho com vencimento em JAN03, onde o preco 
passou, em 120 dias, de R$ 16,20 a R$ 27,00. Se uma operacao fosse realizada nesse 
periodo, para o contrato de venda, os custos associados a encargos financeiros chegariam 
a 9,85% do valor inicial do contrato, para um agente descapitalizado, e 3,12% para um 
agente capitalizado (Tabelas 5.3 e 5.4, pag. 60). O IRPF pago seria nulo nesse caso e os 
custos operacionais totais seriam de 4,18% e 12,2%, para os agentes capitalizados e 
descapitalizados, respectivamente. Nessa mesma situacao, para o contrato de compra, 
teria havido um recebimento de juros equivalente a 3,57%, e um pagamento de IRPF 
equivalente a 15,06%, com um custo operacional total de 13,6% do valor inicial do 
contrato (para todos os agentes). Se os agentes mantivessem suas posicoes pelos 120 
dias, teriam sido esses os custos observados, que curiosamente sao diferentes 
(assimetricos) para posicoes vendida e comprada, especialmente se as operacoes 
envolvem agentes com diferentes graus de capitalizacao e exposicao a tributacao. 

Um aspecto importante observado e a incerteza expressiva que existe no infcio da 
operacao quanto ao valor dos custos ao seu final (Figura 5.9). Os resultados mostraram 
que os custos sao incertos no infcio da operacao e podem crescer quando aumenta o 
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numero de dias da operacao e a volatilidade. Dependem tambem da natureza do processo 
estatfstico que caracteriza a varia£ao de precos (para uma mesma volatilidade e numero 
de dias). Essa incerteza aumenta com volatilidades maiores e operacoes mais longas, se 
traduzindo num risco de custo , um conceito aparentemente novo e particularmente 
importante para as operacoes no Brasil. 

O risco de custo pode, ao lado do risco de base e outros custos de transa§ao, dificultar 
operacoes pelo impedimento de um hedge que possa efetivamente fixar o preco (e 
margens em operacoes). Esse problema seria menos serio se custo fosse fixo, conhecido 
exatamente no infcio da operacao, algo que em geral nao ocorre. 

Ha algumas implica§oes economicas importantes decorrentes de custos expressivos, 
incertos e assimetricos entre agentes, em operacoes em mercados futuros. Custos 
elevados e incertos tendem a reduzir o beneffcio de uma operacao de fixa§ao de precos 
via mercados futuros, inibindo negocia§oes e favorecendo a fixa§ao de precos por outros 
mecanismos que podem garantir precos (e eventualmente adiantar recursos) com certeza 
e menos tributados: CPRs, fixa§ao de precos junto a traders e processadores, etc. Os 
precos fixados por esses outros mecanismos, contudo, incorporam desagios que podem 
ser, em parte, motivados pelos custos e incertezas existentes na fixacao de precos pelo 
mercado futuro. Agentes economicos mais capitalizados e/ou menos expostos a tributes, 
podem ser mais estimulados a utilizar o mercado futuro que agentes descapitalizados e/ou 
mais expostos a tributa§ao. Essa situa§ao e de uma certa forma perversa, pois mostra que 
ha barreiras de acesso ao mercado futuro (via custos de opera§ao mais elevados) para 
muitos agentes que potencialmente poderiam se beneficiar bastante dessas opera§6es (ex. 
produtores rurais do tipo pessoa ffsica). Para trabalhos que estimam a razao otima de 
hedge as implica§6es podem ser grandes, em decorrencia de que em geral assumem 
custos baixos e fixos nas opera§6es. Um entendimento mais completo das implica§6es 
economicas desses resultados, contudo, exigira um esfor§o de pesquisa que foge ao 
escopo deste trabalho. 
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Deve-se destacar que na literatura internacional os custos operacionais em mercados 
futuros tendem a ser praticamente desconsiderados em funcao de sua pequena magnitude, 
quando comparados a outros custos (veja, por exemplo, Nelson, 1985). Essa 
pressuposicao nos parece ser pouco realista para o caso brasileiro. Cabe discutir aqui 
algumas razoes que motivam esses custos elevados e incertos, em magnitude muito 
superior a observada em outros pafses. 

A primeira razao deve-se a: (a) taxas de juros elevadas praticadas no Brasil, em geral 
muito acima das taxas praticadas em outros pafses; e, (b) diferenca significativa existente 
entre as taxas para aplicacao e captacao. Esses dois fatos encarecem as operacoes 
significativamente, especialmente para os agentes economicos descapitalizados, que 
precisam captar recursos para manutencao de suas operacoes de fixacao de precos. As 
taxas de juros afetam os componentes de custo operacional associado a margem de 
garantia, gestao dos ajustes diarios e do fluxo de caixa em geral. Os custos associados a 
margem podem, contudo, ser conhecidos com muita certeza no infcio de uma operacao, o 
que nao acontece com os custos decorrentes de encargos financeiros dependentes do 
fluxo de caixa, que sao incertos. Numa posicao de venda, uma trajetoria de precos 
ascendente pode motivar uma substancial demanda por recursos para fazer frente a 
ajustes necessarios. Em situacao analisada neste trabalho, correspondente ao contrato do 
milho com vencimento em JAN03, onde o preco passou, em 120 dias, de R$ 16,20 a R$ 
27,00, os custos associados a encargos financeiros chegariam a 9,85% do valor inicial do 
contrato, para um agente descapitalizado, e 3,12%, para um agente capitalizado, para 
manutencao do hedge (Tabelas 5.3 e 5.4). As distorcoes existentes nas taxas de juros 
acabam tendo como a conseqiiencia perversa o afastamento do mercado futuro dos 
agentes economicos mais descapitalizados, que em muitos casos sao os que mais podem 
se beneficiar das operacoes nesses mercados. 

A segunda razao deve-se a tributacao. No caso do imposto de renda, analisado neste 
trabalho, em operacoes envolvendo pessoas ffsicas, o tributo e elevado e produz 
distorcoes. Da forma que e definido, o imposto acaba motivando uma parte consideravel 
do custo operacional e da incerteza em uma operacao, favorecendo com isso a utilizacao 
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de outras formas de fixacao de precos que nao sao tao severamente tributadas, e mesmo a 
sonegacao. Como o imposto depende da trajetoria seguida pelos precos, contribui 
significativamente para o componente incerto do custo da operacao, aumentando a 
incerteza com relacao ao custo operacional observado. Nao e nada claro o beneficio 
economico existente para justificar a tributacao de operacoes atraves dos mecanismos 
existentes, que acaba inibindo operacoes que sao importantes para reduzir riscos em 
muitas situacoes. 

Os procedimentos desenvolvidos podem ser estendidos para a analise de outras situacoes 
envolvendo agentes que sao pessoa juridica, com diferentes regimes de tributacao. Os 
custos operacionais para esses agentes podem ser diferentes em funcao de que os 
mecanismos de tributacao utilizados em cada caso tern caracteristicas muito peculiares. 
A extensao dos procedimentos para esses casos e outros casos, contudo, e assunto para 
futuros trabalhos. 

Finalmente, um merito importante deste trabalho e sugerir, atraves dos procedimentos 
propostos, que o custo de operacoes no mercado futuro, no Brasil, pode ser muito 
superior a valores praticos frequentemente discutidos (ex. 0,64%), que tendem a 
considerar simplesmente as taxas da bolsa (TOB, registro e emolumentos). Alem da 
pesquisa apresentar uma estimativa desses custos operacionais, argumentando que podem 
ser elevados, mostra tambem que sao incertos no infcio da operacao, caracterizando um 
risco de custo que juntamente como o risco de base podem dificultar operacoes em 
mercados futuros. Ao considerarmos que as 2 contrapartes da operacao estao sendo 
oneradas, o custo transacional associado a operacao (podem existir outros) acaba sendo a 
soma dos custos operacionais de cada contraparte, o que amplifica o problema. Um outro 
aspecto interessante e que para uma dada trajetoria de pre§os os custos operacionais 
podem ser assimetricos para a posi§ao vendida e a posi§ao comprada. Esforcos na dire§ao 
da mitigacao desses custos e de sua incerteza podem favorecer um maior interesse pela 
utiliza§ao dos mercados futuros como mecanismo de fixa§ao de precos futuros. 
Iniciativas como o FIF e o uso do Funcafe para financiamento de margens e ajustes (no 
caso do cafe) podem ser importantes para a reducao do custo operacional. Com rela§ao a 
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tributaijao pelo IRPF, e diffcil entender a justificativa economica para o mecanismo 
utilizado, que pode introduzir distorcoes economicas desnecessarias e inibir operacoes 
nos mercados futuros importantes para a reducao de riscos. 
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Figura A.l Algumas trajetorias simuladas a partir do modelo estimado para MAI02 
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Figura A. 3 Algumas trajetorias simuladas a partir do modelo estimado para MAR03 
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Figura A.4 Algumas trajetorias simuladas a partir do modelo estimado para MAI03 
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Implicates de custos elevados e incertos para a razao otima de hedge 

Dentro de uma operacao de hedge de venda, o resultado obtido e definido por 

R(X) = P t s + X(P 0 f - Pf ) - CT(X) (4.21) 

S f 

onde P t e P t sao os pre§os observados no mercado fisico e no mercado futuro no 

penodo t, Pq e o preco do contrato no imcio da operacao, X 6 a razao de hedge e 

CT{X) sao os custos de transa§ao (aqui representados pelos custos operacionais), que 

sao dependentes de X. Dentro dos desenvolvimentos realizados a seguir, o 
componente de custo sera representado por 

CT(X) =X(C +ht )c ) + C f I^ 0 (X) (4.22) 

onde c representa a esperan§a matematica do custo e ^ c um desvio aleatorio, com 
esperan5a matematica zero e variancia 1, h e uma constante, l^o(X) e uma fun§ao 
indicadora que assume valor 1 quando X ^ 0 e zero quando X = 0 , e C/ e um 
componente de custos fixo dependente na realiza§ao de algum hedge. A variancia de 
CT(X) e definida por X 2 h 2 . Essa formula§ao tern o objetivo de decompor o custo 
em componentes fixos e incertos, sem perda de generalidade, de forma que os 
impactos desses 2 componentes na razao otima de hedge possam ser mais facilmente 
identificados. 

Para um agente cujas prefercncias podem ser modeladas pelo paradigma da esperan§a 

da utilidade, a razao otima de hedge, representada aqui por X , 6 o valor (possfvel) 
de X que leva o valor de EU(R(X)) ao maximo, ou seja, 
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A* = max 1 EU (R(k)) (4.23) 

onde EU(R(k)) e a funcao que representa as preferencias do agente e &x e o 
conjunto que define os valores possfveis de A . 

Caso 1: Hedge sem custos e sem risco de base 

A situacao mais simples possfvel com relacao a razao otima de hedge e caracterizada 
pelo caso em que os custos da operacao sao nulos e inexistem riscos de base, ou seja, 

o preco recebido na transa§ao no mercado fisico, P t , e identico ao pre50 recebido no 
mercado futuro e P t 



Se assumirmos que as preferencias dependem somente da esperan§a e da variancia do 
resultado, ou seja, 

EU(R(X)) = E(R(X)) - ^V(R(X)) (4 . 2 4) 

onde r e um coeficiente de aversao ao risco absoluto (do tipo Arrow-Pratt). Os 
paragrafos seguintes exploram 3 situa§6es distintas quanta a natureza dos custos 
operacionais. Na primeira situa5ao, assume-se que o custo e nulo. Na segunda 
situacao assume-se que o custo e fixo, igual a c (nesse caso h=0). Na terceira 
situa5ao assume-se o modelo completo, com custos nao nulos e incertos. Para facilitar 

f f 

a exposi§ao, sera assumido tambem c / = 0 (custos fixos nulos) e E(P t ) = P 0 , ou 

seja, o pre50 no periodo 0, no mercado futuro, e um estimador nao tendencioso do 
pre§o no periodo t. 

Razao otima de hedge com custos nulos 

Com custos nulos temos uma simplifica§ao de (4.21) por 
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R(l)=P;+MP 0 f -P/), (4.25) 

e preferencias representadas por (4.24), chega-se a 

EU(R(X)) = E{P} ) - ^(c 2 + XW f - 2Xa sf ) 

A razao otima de hedge sera o valor de X que soluciona (4.23). 
ser obtida pela analise das condicoes de primeira ordem, dada por 

Xa 2 -a sf =0 

As condicoes de segunda ordem, obtidas a partir da derivada segunda de (4.26), com 
relacao a X, que e estritamente negativa nesse caso, garantem ser (4.26) uma funcao 
estritamente concava em relacao a X . Com isso, qualquer maximo local sera tambem 
um maximo global, atendidas as restricoes a respeito dos valores possrveis para X, . O 
resultado das condicoes de primeira ordem, (4.26) 

EU{R(k)) = E(P s f )-\c--(<5 2 s + X 2 a 2 + X 2 h 2 - 2Xa sf - 2Xha s ^ + 2X 2 ha^) 



(4.26) 

Essa solucao pode 



o mercado futuro preferencias do agente com relacao ao hedge poderao ser 
expressas, a partir de (4.21) e (4.22), por 

EU(R(X)) = E(Pf)-Xc-^(G 2 +X 2 a 2 +X 2 h 2 -2Xa sf -2Xha s ^ + 2X 2 ha^) {42 5) 

onde cr^, a 2 e Jih 2 representam, respectivamente, as variancias do preco no 
mercado ffsico, futuro e de CT(k) , e G s f,O s ^ e a ^ representam, respectivamente, as 
covariancias duas a duas do preco no mercado ffsico, futuro e de ^ c . 

Os paragrafos seguintes 
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Atraves das condi§6es de primeira ordem, a razao otima de hedge sera obtida por: 

X * = 2 „, <<V +k a 4 

O f +h z + 2h G£ * r 

^ = ^,2 ^ n „ (®sf + h P*?s ~ 7) (4.26) 

<5f +n +2n P ftp f r 

caso , ou um dos pontos extremos de &x , caso A, e £\ . Em particular, 

quando £\ = [0,1] , os pontos extremos correspondem as situacoes em que nao se faz 
hedge (k = 0) ou se faz o hedge completo (A- = 1) . 



Se os custos sao fixos (sem incerteza), a expressao (4.26) poderia ser simplificada 
para 

* _ O s f _ G s 
-~2"-P#7~ (4.27) 

que e a solucao teorica para o "a razao de hedge que leva a variancia minima", na 
ausencia de custos de transacao. 
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